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ABSTRACT 


Starting from Breit’s relativistic equation for a system of two electrons, 
it is shown that for a hydrogen molecule (or for a system of two electrons 
moving in a field of cylindrical symmetry) the component of the total 
angular momentum (J,) along the axis of the molecule (axis of symmetry) 
is a constant of motion. Thus every eigenstate of the system is simul- 
taneously an eigenstate of J, also, and a state of the system will specify, 
besides its energy, only the eigenvalue of the component of the angular 
momentum parallel to the axis of symmetry. The form of the four large 
components of the wave function relating to their dependence on the 
azimuthal co-ordinates has been given. 


The case of Russel-Saunders approximation has been considered 
in detail and the nature of the components of the wave function for the 
singlet and triplet states has been discussed. It is shown that the wave 
function for the ground state of the hydrogen molecule could be expressed 
as a sum of a set of symmetric functions of which the first term is the 
Heitler-London function, and that the wave function for a triplet state 
should be a superposition of anti-symmetric molecular orbitals. It is 
shown that relativistic theory brings about in a natural manner the facts 
relating to the ground state of the molecules C, and O,. Finally, some 
remarks are made concerning the case of molecules for which the spin- 
orbit and the spin-spin couplings are strong. 


I. INTRODUCTION 


THE theoretical solution of problems of molecular chemistry are based usually 
on one or the other of two methods, known as the Molecular Orbital Theory 
and the Valence Bond Theory.! In the former method, the binding electron 
is supposed to belong to the entire molecule and is assumed to be moving in 
the average potential field of the other electrons and the nuclei. The latter 
is the well-known theory of hydrogen molecule given by Heitler and London. 
An important consequence of this theory is that the ground state of the hydro- 
gen molecule is the singlet state in which the spins of the two electrons are 
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opposite to each other and that the states with the spins parallel to each 
other are unstable. In both the above-mentioned methods, the spin of the 
electron is introduced as an independent extraneous factor in such a way as 
to comply with the requirements of the Pauli Exclusion Principle, but is not 
brought out as an integral part of the theory itself. In view of the very 
important role that the spin of the electron plays in the chemical binding in 
molecules, it is desirable to develop a theory based on a relativistic equation 
which will bring out the spin as an essential part of the mathematical for- 
malism itself, and to analyse the consequences of such a theory. This, pre- 
cisely, is the object of the present paper. 


Il. THE HAMILTONIAN OF A SYSTEM OF TWO ELECTRONS 


Let us consider a system of two electrons interacting with each other. 
Let ¢', A’ and ¢", A™ denote respectively the scalar and vector potentials 
of electrons I and II disregarding their interaction. A relativistic wave 
equation for a system of two electrons interacting with each other and with 
an external electromagnetic field was given by Breit? in the year 1929. 


The equation is 


Pot J (ax'pr’ + ax" pel) + (aq! + ay!) me 


k=1,2,3 
e tay 4 (a! +B) (alt. r) ? 
+ Yer | *%' 8 ee ee ]. 
where 
pore” é I. 
Pk — ih Xk! + (¢) Ak > 
ie ae a ! e II 
Pel = — ih oat + (<) Ax 
and 
th , ec ON 2 
Po “e+ @ + $") = (2) 


The matrices a;,' and «;!! are the usual Dirac matrices, but now they operate 
on the first and second suffixes respectively. The representations that we 
choose for these are given in Appendix I. The wave function % has sixteen 
components which will be designated by %mn (m, n = 1, 2, 3, 4), the first 
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suffix m referring to electron I and the second suffix n referring to electron II. 
If a! and 5" be any of the above-mentioned matrices, then 


(2'D"P)m,n = E Amx'bnr fe (3) 


Equation (1) is thus equivalent to sixteen differential equations in the sixteen 
functions %mn. 


By writing in (1) 

@, = P1%1;5 As = P1%; ag = P1%3 (4) 
where p, 0,, 02, 03 are the matrices (4) of Appendix I, we see that the 
Hamiltonian of the system is given by 

H = e(¢' + $1) +% —pilc(!. p') 


cs pric (a! . p") ‘eit mc2 (a,! + a4!) 


2 1 1 
ie 55 | (ove » pyle!) + (p;'o' . r) (pre es =] (5) 


r? 


Ill. THE x-COMPONENT OF THE TOTAL ANGULAR MOMENTUM 


We shall consider the case of the hydrogen molecule which is the simplest 
example of a two-electron molecular system. Let the origin be chosen at 
one of the nuclei and let the line of action of the molecule coincide with the 
positive direction of the x-axis. If the nuclei are denoted by the suffixes A 
and B, the potentials ¢' and ¢"' are given by 


! 1 
I — ae 
$ =e(-- +5, 


and 


poet re 


Ta, lb, 


where ra= Vx*+y?+2* and ro= V(x—R+y*+ 2%, R being 
the internuclear distance. 


Let m,' and m,| denote the x-component of the orbital.. angular 
momenta of electrons I and II respectively. Then the dynamical variable 





4 K. S. VISWANATHAN 


Mz = (mz' + m;") represents the x-component of the total orbital angular 
momentum of the system. In the same way, the quantity 





Sx —_ 5 (ox! + ox") 


gives the component of the total spin angular momentum parallel to the axis 
of the molecule. The sum J. of these two (i.e.) 











Je = mz! + ma! +5 (ox! + oz") 6) 
represents therefore the x-component of the fofa/ angular momentum of the 
system parallel to the axis of the molecule. 

It has been proved in Appendix II that 
J<H — HJ; = 0 





(7) 


One can see from the Appendix that neither S; nor M, commutes with the 
Hamiltonian but only their sum has the property of commuting with H. 





Thus the component of the total angular momentum parallel to the axis 
of the molecule is a constant of motion. 





Since J, commutes with H, every eigenfunction of H will be simulta- 
neously an eigenfunction of J, also. An eigenstate, thus, is specified by two 
parameters, namely, the energy of the state and the eigenvalue of the com- 
ponent of the total angular momentum parallel to the axis of the molecule. 








IV. THe EIGEN #’s OF Jz 


Let us denote the eigenvalues of J,, by M&# where M is yet undetermined. 
The equation 


J = Mik 








(8 a) 






[ ma + m," + 5 (on! + oz") | ob = Mad (8 5) 


will then be equivalent to sixteen equations in the sixteen components of 4. 
We write here only four of these equations, namely, those involving the four 
large components %g, Yg4, Yyg and 44, of #. From (3), it follows that the 
(3, 3), (3, 4), (4, 3) and (4, 4)th components of the equation (8) are given by 
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— 1 (56; + ag) #9 = OM — Dt 


— (3p, + agg) Ye = Mil 
— (55, + 3g.) #0 = Mito 


Ee (sr + x) hag = (M + 1) day 


Yas = Soa (115 Fo, 91, 99) ef (m,-19,+m,9,) 

Psa = Sa (as Tay 9, 99) ef M4492) 

has = Sas Tr» Fay 91, 99) eb 9.4.92) 

toss = faa Wry Fey 94, Og) ef 9:4 419,) (10) 
where m, + m, = M. 
Since the ’s should be single-valued, the numbers m, and mg, are integers. 


One can verify that the other twelve small components satisfy differential 
equations of the type (9) and their form is therefore similar to the expressions 
given by (10). 


Now it has been shown by Breit?: * that equation (1) involving sixteen 
components can be reduced to an approximate equation with four compo- 
nents expressed in terms of Pauli’s spin matrices. This reduced equation is 
given by (6) of Appendix I and is equivalent to a set of four equations in the 
four large components of y% which, following Breit, we denote by 3, 34, 
ag and 4, respectively. Since the equation (6) of Appendix I is sufficient 
for all physical purposes and gives quantitative results correct to the order 
of (v/c)?, we shall in the following sections be discussing the consequences 
following from this only. 


The matrices o! and o!! in (6) of Appendix I are the Pauli matrices and 
they operate on the suffixes 3 and 4. The operation rule for the product 
of any two matrices on ¢ is still governed by (3), but now the summation for 
k or / extends over the numbers 3 and 4 only. One can verify that the 
equation J, = Mim is identical with the four equations given by (9), and 
the form of the four large components is therefore still given by the 
equations (10). 
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By solving the simultaneous differential equations (7), Appendix I, one 
can determine the functions 55, 54, %43 and #44, respectively. Thus a state 
ws which we denote by the column vector [is, 34, 43, 44] Consists in general 
of four non-zero components. 


V. THe Russet-SAUNDERS COUPLING 


When the spin-orbit interaction is small compared to the electrostatic 
interaction of the two electrons, the terms involving B', B", Cz as well as the 
last four terms in the expression for A in equation (6, Appendix I) could be 
ignored in the first approximation. This is the assumption underlying the 
case of the well-known Russel-Saunders coupling. The equations (7, 
Appendix I) now reduce to 


(H — W) #53 = 0; (H — W) #54 = 0 
(H — W) 443 = 0; (H — W) $4, = 0 (11) 
where H is the non-relativistic Hamiltonian of the system. The above equa- 


tions are satisfied by four different states of the system in each of which all 
components, excepting one, are zero. Thus, solutions of (11) exist which have 


the form 
$y = [yg3, 0 0, 0]; fo’ = [0, #54, 0, 0] 
bs’ = [0, 0, 4s, 0]; %, = [0, 0, 0, bya] 


Since the spin-orbit coupling terms are neglected in the Hamil- 
tonian (6, Appendix I), one can see that both 


Mz = m,' + mz" and Sz = ; (o,' + oz") 


are constants of motion in this case, and further 4, and 4, are states belong- 
ing to the eigenvalues # and — hi of Sz, respectively, whereas the states %,’ 
and 7,’ have zero as the eigenvalue of S;. From the form of #5, and 4; 
in (10), one can see that the states 

fo = pe’ + ty = [0, vga, Yas, 0] 
and 


bs = Pa — ty’ = [0, tea, — Yaa, 0] 


have the eigenvalue M for the x-component of the total orbital angular mo- 


mentum and zero for the x-component of the total spin angular momentum. 
Fucther if 


>? — 4 [(c,! + ox!')? aa (cy! + oy")? + (cz! + oz'!)?] 
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one can verify from (3) that the states 4,, 4. and #, are eigenstates of =? 
belonging to the eigenvalue 2 (or the value | for 2), and the state ws, has the 
eigenvalue zero for 2. The former set of states are thus the well-known 


triplet states of two electron systems whereas the state 7, represents the singlet 
state. 


Now according to the Pauli Exclusion Principle, every state should be 
anti-symmetric with respect to an exchange of the two electrons. An exchange 
operator, which we denote by P, interchanges the spin states of the two 
electrons as well as their positions and thus under its influence, a component 
of the state %mp (1, 2) gets transformed into %nm (2,1). The state 


} = [ibys (1, 2), 34 (1, 2), Has (1, 2), as (1, 2)] 


therefore changes into & = [vgs (2, 1), tag (2, 1), Haq (2, 1), aq (2, 1)] under the 
influence of P. In order to satisfy the Exclusion Principle, the following 
relations should therefore be obeyed by the components of a state: 
(1) 3 (2, 1) = — pag (1, 2); (2) Haq (2, 1) = — pas (1, 2) and (3) Yaa (2, 1) 
= — %,4,(1, 2). Onecan see that the states , and 4, will be anti-symmetric 
if the components ¢,; and y44 are anti-symmetric functions of the position 
of the two electrons. Since #4, and x4, are solutions of the same 
Schroedinger equation in the Russel-Saunders approximation, the state », 
will be anti-symmetric if its components 3, and 43 are such that #4 (1, 2) 
= hyg (1, 2) = — oq (2, 1) = — ag (2, 1), and similarly the state 4, will be 
anti-symmetric if 


tsa (1, 2) = Ysa (2, 1) = — as (1, 2) = bas (2, 1) 


(a) The symmetric and anti-symmetric  states—Let nim (Ta,) 

denote a complete set of eigenfunctions of a hydrogen atom whose nucleus 
is at position A. Then these functions form a complete basis for the expan- 
sion of any arbitrary single-valued function of the position r, = (4, y1, 21) 
of the electron 1. Similarly, the set of eigenfunctions (%yym (fp,) of a 
hydrogen atom with its nucleus in position B form a complete set of eigen- 
functions for the expansion of an arbitrary function of the position r, of the 
electron 2. Then, by a well-known theorem,‘ the set of functions 


nim (Ta,) bn'vm’ (rd,) ( 


form a complete basis for the expansion of an arbitrary function % (¥,, r2) 
of the two electrons, that vanishes when r,;—>co and r,—>oco. Thus we 
can write 


 (¥,, ¥2) = 2 Fan’ ¢n (Fa,) bn’ o,) (12) 
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where we have written n and n’ for the set of three quantum numbers (n, /, m) 
and (n’,/',m’) of the hydrogen atoms. 


Now if the wave function ¢ (r,, r.) is symmetric with respect to an inter- 


change of the two electrons, we get from (12) by interchanging the suffixes 
(1) and (2) that 


# (¥,, 2) = 2 Fan’ n (fa,) Sn’ d,) (13) 


The above equation alternatively follows from the expansion of 
# (¥,, f2) in terms of the set of functions %p (Tg,) 4 (Fp,) .... Which also form 
a complete basis for the expansion of a function of the two electrons. 


By adding (12) and (13)* we see that any symmetric function of the 
electrons 1 and 2 can be written as 


p* (F,, Fo) = a FSnn’ [bn (fa,) dn’ (Po,) + Yn Ca,) on’ Fo.) (14) 


where fnn’ =4F nn’. The suffix + has been used to indicate that the func- 
tion % is symmetric with respect to an interchange of the electrons. 


In the same way we can see that an anti-symmetric function J (r,, r.) of 
the electrons has the form 


be (71, To) = 2 Lnn’ [Yn (Fa,) dn (Fo,) — bn (a,) dn’ (rp,)] (15) 


(b) The singlet state——We have seen before that the singlet state x, is 
represented by the vector [0, %,;*, — 54°, 0] where .,* is a symmetric func- 


tion of the two electrons. Now from the equations (10) it follows that 5, 
has the form given by 


Ysa = Soa 11> Ta, 91, 92) et (M,9,+M,9,) (10 a) 
where m, and m, are integers satisfying the relation m,-+m,=M. Thus 
the general form of #4 consists of a superposition of terms of the type (10 a) 
for all integral values of m, and m, satisfying the condition m, + m. = M. 

Let us consider the case where the component of the total angular 
momentum parallel to the axis of the molecule is zero, i.e., m, + m,=0 


or m, = — m,. It is well known that the ground state of the hydrogen mole- 
cule belongs to this species. Thus from (14) and (10 a) it follows that the 
eigenfunction of the ground state is given by 


Yea (W182) = LF F fan [enim (Fa,) dn'v-m (Fo,) 


nim a’ 


+ tnim a,) dnv—m o,)] (16) 


* We assume the absolute convergence of these series. 
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The first term of the above series is 


[¥6100 (Fa,) 100 (Fb,) + +100 (Fa.) ¥4100 (,)] (17) 


which is the well-known Heitler-London wave function. This function 
explains qualitatively in a satisfactory manner the chemical binding of the 
hydrogen molecule, but quantitatively it gives a value of 3-14e.v. for the 
binding energy of the molecule, the observed value being 4-72¢.v. The 
quantitative inaccuracy of the Heitler-London function is not surprising 
when we recall that it is only the first term in the expansion of the wave 
function in terms of symmetric wave functions built from the hydrogen atomic 
orbitals. Since the discrepancy between the energy value derived from (17) 
and the experimentally observed one is not very large, one can hope that 
the series (16) converges rather rapidly and that the addition of a few more 
terms to (17) will lead one to the correct value for the binding energy of the 
molecule. 


(c) The triplet states—As an example of the triplet states, consider the 
state %, which is described by a non-vanishing anti-symmetric component 
sg (¥;,%2) and three vanishing components. From (10) we have 


bsg = Sag (11 Te, 91, 99) et (M,$,4+M,9,) (18) 


where m, + m, = (M — 1). 


Since the component of the total spin angular momentum for this state 
is A, the component of the total orbital angular momentum is equal to 
(M — 1)%. The general form of 4, will consist of a superposition of terms 
of the type (18) for all integral values of m, and mz, satisfying the relation 
m, + m,=(M — 1). Thus from (15) it follows that 4,,- is given by 

sg (KF, To) = ~ = Snn’ [Ynim (a,) Yrt'm’ (Fd,) 
— nim fa,) bn'v'm’ (fd,) (19) 
where m +- m’ = (M — 1). 


In other words, the state 4, is a superposition of all anti-symmetric mole- 
cular states for which the component of the total orbital angular momentum 
is equal to (M — 1)’. In the same way we can obtain the wave function 
for the state 4, by superposing anti-symmetric molecular states have the 
eigenvalue (M+ 1) for the component of the total orbital angular 
momentum parallel to the axis of the molecule. 


Regarding the state %, we shall consider the particular case M = 0. 
ws, is described by the vector [0, 45,,%34¢,0] where %,;> is anti-symmetric 
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with respect to an interchange of the two electrons. From (15) and (10a) Ji 
it follows that the form of ,;- is given by int 
Yq (1, Fo) = = - Snn’ [Lnim Ca,) &n'v—m (Fo,) . 

— tntm (a,) dart (Fo) Co ew 

The first term of the above series is of 
[100 (Fa,) Ya00 (OF,) — Yro0 CaF.) P00 (F,)] (21) - 

It is a well-known fact that this function leads to an unstable chemical state. sl 
fc 


(d) The case of C, and O,.—From (19) or (20), it follows that the wave 
function for a triplet state is a superposition of a number of terms each of 


whicn is anti-symmetric with respect to an interchange of the two electrons. ¢" 
Further, all the terms in a series belong to the same species, i.e., the states tl 
belonging to these terms have all the same eigenvalue for the x-component A 
of the total orbital angular momentum. C 


According to the Heitler-London theory, the ground state of the hydro- 
gen molecule is the singlet state and the triplet states are unstable. The 
theory therefore requires as an essential condition for a stable bond that the 
spins of the two electrons should be antiparallel. This rule is obeyed by a 
large class of molecules, but nevertheless, exceptions are also present of which 
we may cite the exampks of C, and O,. The ground state of C, is the state | 
designated in spectroscopic notation by *J7, whereas for oxygen, the ground 
state is *2-. Oxygen is in addition paramagnetic. Heitler® and Nord- 
heim-Poschl’ showed that by adding to the wave function (21) a few more 
terms of the same species, or in other words, by taking into account the inter- 
action between states of the same species, the deviations between theory and 
experiment can be removed. Now in the Russel-Saunders approximation, 
= has the value 1 for the triplet states and 2, has the values 1, 0, — 1 res- 
pectively for the states %,, 4. and %,. Further each term in (19) or (20) cor- 
responds to the same value of the component of the total orbital angular 
momentum parallel to the axis of the molecule. The expansion (19) is, there- 
fore, nothing but a superposition of states of the same species. Relativistic 
theory offers thus a clear explanation of the facts relating to C, and O, and 
brings about naturally the concept of interaction of states of the same species. 


VI. StrRonG Sprn-OrsBiTt COUPLING 


The discussion of the last section holds good only for the case of Russel- 
Saunders coupling. Strictly speaking, the Russel-Saunders coupling and the 
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j-j coupling are two extreme ideal cases and the coupling in most cases is 
intermediate. When the spin-orbit interaction is strong compared to the 
electrostatic interaction of the two electrons, the above discussion will break 
down. In such a case, only the component of the total angular momentum 
along the axis is a constant of motion and none of the quantities 2, 2, and 
M., is a constant of motion. An eigenstate of the energy will be an eigenstate 
of J,, also and conversely every state of the system will specify, besides the 
energy, only the x-component of the total angular momentum. For the case 
of strong spin-orbit coupling, the usual spectroscopic notation which pre- 
sumes that all the quantities 2, 2, and My, are constants of motion, is there- 
fore inconsistent with theoretical considerations and requires. modification. 


When the spin-orbit coupling is strong, the solution of the equations 
(7, Appendix I) will lead to non-zero values for all the four components of 
the wave function. Further it follows from equations (6) and (7) of 
Appendix I that a state will be anti-symmetric with respect to an interchange 
of electrons only if the components have the following properties: 


(1) 33 (2, 1) = — og (1, 2); 
(2) gq (2, 1) = — pa (1, 2); 
(3) tag (2, 1) = — Haya (1, 2). 


Thus 3 and 44 are anti-symmetric functions, while 4, or %4, can be a mix- 
ture of both symmetric and anti-symmetric functions. The form of the 
former components is therefore given by (15) while the latter can be expanded 
as a series of the type (12). Now the energy of the system is given by 


_ J b*Hdbdr,dr. 


W =" Pender, = 


where s* is the transpose of the vector % and is given by the row vector 


[bo3*, ba*, tas *, taa*]. 


In practice, the computation of the energy will be facilitated if (1) the 
series for the function #33, etc., are terminated at a particular point and their 
coefficients determined by the variational principle and (2) if use is made 
of the tables giving the values of exchange, coulomb and other related 
integrals.*®: 9 1° 


We can now interpret the quantities 4,7, Wg4? + y43? and 4,4,” as the 
probabilities of finding the spins of the electrons respectively (1) parallel to 
the x-axis, (2) opposite to each other and (3) parallel to the negative direction 
of the x-axis. The ground state of the molecule will be the state having the 
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least value for the energy of all the states, subject to the condition that it is 


The author’s grateful thanks are due to Professor Sir C. V. Raman 
for his kind interest in this work and valuable suggestions. 
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The Relativistic Theory of Chemical Binding 
APPENDIX I 


For the sake of easy reference, we collect here some of the results and 
equations of Breit that are made use of in the paper. The relativistic wave 
equation of Breit is given by 


Pot Z (ax'pe! + an" pelt) + (ay! + 24!) me 
ke1,2,3 

~=0 (1) 

e” toy. tt 4. (a! + ¥) (alt - F)) 
+ dap {Z ae! += 7? j 
Here pp, p' and p" are given by equations (2) of Section II. The wave func- 
tion % has sixteen components which are denoted by %mn (m, n = 1, 2, 3, 4). 
The matrices a! and o" are the usual Dirac matrices, but now they operate 


on the suffixes I and II respectively. They are given by 


0 
0 
| 
0 
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0 —-i 0 0 
ee he 2 gage ee 
ee 9 a \ 


0 
0 
1 
0 


- ©o& Oo O&O 


eo ¢o8 @ = 


ed 


equation (1) becomes 


(» + pi! (o!. p) + py (oe! . p!) + (a4! + a4!) me 
» pa!) + (p;'a' . r) (pyital' - 1) Ate ©) 
r2 


\ 


e2 
+ 2cr [ v0 


It follows that the Hamiltonian of the system is given by equation (5) 
of Section II. 


Next the Breit equation in the Schroedinger-Pauli form is given by 


JE — 2mc? + eV — _ ((p)2 + (p")2] + (p')* + (p")* 


8m3c3 


e2 
TF apace LO (pts BE) 178 Gait — x4!) Gait — 255") pil 5" ] 


sb <=. ((H', a’) J. (H", a)] + =e ((é', p) + (64, p")] 
p' p" 
+ eam [(e x m’ ot) + (en 4 at) 
1 yl i — I 
+ 2e(*— mn xP ie t) + 26 — ec, 
2 
~ (452) (V7 "o") (V'o") ryt p= (6) 
where 

H! = curl A', H"™ = curl A" 
é! = — gradi V = — (3. = yh st) V3 


I I e 
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() means the differential operators p', p" apply only within the { ). 
The o’s are the Pauli matrices and have the representations given by 


an(h Sia en? 4) 


Let us write 
E — 2mc? = W; 


(p')* + (p")8 


l 2 
A = — 5 (P+ (WIS + ras 


+ act iP p . p" + Sr-3 (x! — x;') (x;" = x;') pi’ p;"] 
, Chi 
Sam2c2 Ce BY) + C6", p's. 
eh p' ri — rll pil 
B= Samc? [ e x me de 7 a = (Bz, By’, Be’) 
eh p! rl! — ri p! 
BU — Same? [ 6" x = + 2e _— x F| = (@,", B,"5 Bz") s 
and 
eh \? 
— (45) (0 DT 08) ry = 2 Cayos'oy!" 
where 


Cry =-—(455) = In! dy ve 


Then the equation (6) is equivalent to. the following four equations: 


(A + By! + By! + Cox + W) ts + (By —iBz" + Coy -- iCzz) ta, 
+ (By’ =n iBz' + Cyx — iCzz) bas + (Cyy = Czz 
oa iCyz so Czy) $a, = 0. 
(B," + iB," + Cry + iCyz) #33 + (A + Bz By" — Cox + W) tsa 
+ (Cyy + Czz) bas + (B,' — iB, — Cyx + iCzz) Baa = 0 
(7) 
(By! + iBz' + Cyx + iCzq) tg + (Cyy + Czz) bea 
+ (A + Bet — Byt — Coz + W) das + (By"— iB,"— Czy 
+ iCzz) Baa =0 
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(Cyy —C2z+i Cyz + Czy) sg + (B,' +i B,* ~ Cyx i iCzx) Psa 
+ (B," + iB," — Cry — iCyz) t43 + (A — B," i B,." 
+ Crx + W) Pas = 0. 


APPENDIX II 


Proof of the relation J,H — HJ, =0.—We have seen that the Hamil- 
tonian for a hydrogen molecule is given by 


2 
H = — e(g + $) + < — pyle (ot . pt) — pyc o™ - p®) 


° 


. e , (po! .r Hl , r)) 
— me? (af + af) — = j pio! - pol) + pyo' - r) (p,"ta! )) 


2 r? 
(1) 
where 


eS and gt=e( +1); 


ra, ro, 


=e 


r is the interelectronic distance. Each of the angular momentum variables 


m,' and m,," commutes with rg,, ’b,, "a, and rp, and hence with the potentials 
¢@ and ¢". Further 


[mg! + my%, 1?) = [mz! + my", Z xy? + 2x? — 22 x49] 
= — 2(ViZ2 — YozZ1) — 2 (Ye21 — ViZ2) = 0 (2) 


(m,' + m,") commutes with r® and hence with any function of r?, particularly 
with 1/r and 1/r°. 


Now 
ih (rite! + my") = (mz! + mz") H — H (mz! + mz") (3) 
Since (m,' + mz") commutes with ¢', ¢", I/r and 1/r®, we have 
ih (tig! + mz") = — cpy' {mz" (a! . p') — (a! - p’) mz) 
— ep," {mg"* (ot « pH) — (ol. pl) mai 


_ + (my! + mz") (p:'o! - £) (p;"4o™ . r) 4 
2r? 1 — (p,fo! . r) (p:Mo™ . r) (my! + al (4) 
Now 


{mz' (a! . pt) — (o . p') mz} = (ol, m,"p' — p'mz') 
= ih (G2p3 — ope)' (5) 








Sir 
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Similarly 
tmz" (a4 . p¥) — (o" . p®) mz"} = ih (ops — o5p2)" (6) 
Next the following Poisson bracket relations can be verified: 
x + mz", (%2 — %) 2 — ma) — 2—%) Ge — %) , 
rs rs ’ 


re r3 


r3 r3 





es 

oe! + mz", (X2 — %) (Z2 - — ~#)] = (x2 — %1) 2-91) : 
~ xz + mz", (Ye — Vi) G— “| “= [‘: — a ~~ — wy" : 
| 


mz! + my", (ye — "| = 2 (Ye — — Ya) (2. - = 2) , 


a r3 


F * + m;,", (2, — 2Y" Peay 2(y2 — Ys) 2s — 2) } (7) 
r 


r3 
Hence if we denote the components of the vectors p,'c! and p,;"o" by 
(az, ay}, az‘) and (az", ay", az"), 
we get from (3), (4), (5), (6) and (7) that 
ih (mhy* + mz") 
= — ih cp;' (o2p3 — o3ps)' — ih cpy" (oops — o3P2)" 


(Az'ay" +- Ay'ay") (X2 — 1) (Z, — 21) 


i — ihe* — (ayz'az" + az'az") (x2 — 1) 2 — Yi) (8) 
2r8 ss (ay'a,™" + az'ay") [(Z2 — 2,)® — (ve — J1)*] 
+2 (ay'ay" — az'az") (Ye — V1) (Ze — 24) 
Next we have 
(cz!o! — —'o,', p') = 2i (csp. — Oe)s)' 
(7gHo! — cog", pl) = 2i (ape — apa)" 0) 


Also 


+ ox", (p,' 7! ° aD — 
r3 
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Further 
og! (pyle? « x) (pH. 2) — (pylo! « x) (pHo™ « nog! 
= 2i {az' (Ye — Yi) — ay’ (Z2 — 21)} (ax! (X2 — Xi) 
+ Gy" (Ye — Wr) + az" (Zz — 23)} 
ox! (pylo! « 1) (pyt¥o « 4) — (po! « 1) (p:!'o% « 2) on" 
= 2i {az" (Ye — Y1) — ay" (Ze — 21)} {ax" (X2 — 1) 
+ ay! (Ye — i) + @z' (Ze — 24)} (11) 
From (9), (10) and (11), we get 
ih (oz' + oz") 
= (o,' + o,")H —H (oz) + oz"! 


= — 2icp,' (csp, — oops)! — icp," (op, — eps) 
(az'az" + az'ay") (v2 — V1) (x2 — x) 
— 2ie? + (az'ay" + az"ay") {(ve — Wi)? — (Ze— 2;)7} 
r3 1 
: + 2 (a@,'a," — ay'ay") (Ye — V1) (22 — 24) 


De (ay'a"" > ay"'az') (Z. — 2,) (x2 — *) 


(12) 


From (8) and (12) it follows that 
ritg! + tite!" +5 (ox! + oz!) = 0. 
Or, in other words, 
Jeg = mz! + mt + # (ox! +o,") is a constant of motion. 


Fields possessing cylindrical symmetry.—Though the above result has 
been proved for the case of a hydrogen molecule only, it holds good in general 
for any system of two electrons moving in a field possessing cylindrical sym- 
metry. If the x-axis is chosen as the axis of symmetry, the potentials ¢1 
and ¢"! in this case are functions of two variables only, i.e., ¢' = ¢' (x,, r,) 
and ¢! = ¢" (x., r.), r,; and r, being the distances of the electrons from the 
origin. Since M, = m,' + m," commutes with all the variables x,, x2, r, and 
ra, it commutes with the first term in the Hamiltonian (1). The remaining 
terms do not involve ¢! or ¢" and are the same as for the hydrogen molecule. 
Thus for a system of two electrons moving in a field of cylindrical symmetry, 
the component of the total angular momentum along the axis of symmetry 
is a constant of motion. Since for most diatomic molecules and individual 
bonds in complex molecules the charge cloud approximately satisfies cylindri- 
cal symmetry, the above result can be applied to such cases also. 
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INTRODUCTION 


Tue structure of L-tyrosine hydrobromide which is isomorphous with the 
hydrochloride has been discussed in the previous paper (Srinivasan, 1959, 
hereafter referred to as Part I). The structure had been solved by the appli- 
cation of the Difference-Patterson (Kartha and Ramachandran, 1955) and 
minimum function (Buerger, 1951) methods. The refinement of the struc- 
ture of the hydrochloride was taken up for the sake of confirmation and 
elucidation of the various features found in the other compound. The 
details are given in this paper. 


EXPERIMENTAL 


The crystallographic and other experimental data are available in the 
previous paper (Part I) and are summarised here for convenience. L-tyrosine 
hydrochloride (C,H,,NO,HCl) belongs to the monoclinic system with space 
group P2,. The unit cell dimensions are a=11-:07A, 5b=9-03A, 
c= 5-09 A, B = 91-8°, V = 508-6 A* and Z= 2. Thecalculated density is 
1-619 g./cm.* while that determined by the method of flotation is | -64g./cm.* 


Photographs for intensity estimation were taken using the multiple film 
technique and the intensities were estimated visually. Aiter correction for 
the L-P and absorption factors, they were placed on an absolute scale by 
Wilson’s method. The coefficient B in the isotropic temperature factor 
exp-(B sin? 6/A?) was found to be 2°35 A? from the Wilson pot. However, 
it was found later that this value had to be altered appreciably, namely, from 
2°35 to3-00 A®. As the cu:pound was rather unstable the needles were 
enclosed in a quartz tube during exposure. However, the data for the [010] 
zone could be collected only at a later stage as it was difficult to enclose the 
crystal cut perpendicular to the needle axis inside the quartz tube. Finally 
the crystallite was dipped in a solution of cellulose acetate to obtain a thin 
film coated on it which prevented decomposition. 
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The atomic scattering factors of Viervoll and Ogrim (1949) were used 
in structure factor calculations. 


The total number of reflexions theoretically accessible using Cu Ke 
radiation was 145 in the [001] zone and 152 in the [010] zone while the num- 
bers experimentally accessible on the photographs were respectively 138 and 


142. Of these 10 in the [001] zone and 12 in the [C10] zone were observed 
to have zero intensity. 


REFINEMENT OF THE b-AxIS PROJECTION 


In order to facilitate the discussions the following notation is used. 
The calculated structure factors usually denoted by F¢ (hk/) is now replaced 
by F‘) (hkl) where the superscript p denotes that the structure factors have 
been calculated for a particular set (corresponding to the p-th stage of refine- 
ment) of co-ordinates listed in Tables I and I]. The same symbol with a 
subscript H denotes that the calculated structure factors include hydrogen 
contributions. The same type of symbols have been used, where necessary, 
for denoting the reliability indices also. Thus, for example, Ry (H0/) 
denotes the reliability index for the set of structure factors Fy,* (h0/) which 
correspond to the 5th set of co-ordinates for this zone and include hydrogen 
contributions. If sets of structure factors are different due to reasons other 
than change in co-ordinate parameters of the atoms (as for example change 
in temperature factor, rescaling, etc.) they are denoted by F,, (A0/), 
F,,%" (h01), etc. The observed structure amplitudes are denoted as usual, 
by Fo (hkl) and the indices hk/ are dropped if there is no ambiguity. 


REFINEMENT OF CHLORINE POSITION 


Since the present compound is isomorphous with the hydrobromide, 
a set of co-ordinates taken from a fairly advanced stage of refinement of 
the hydrobromide would provide us with a good approximation for starting 
the refinement. With regard to the co-ordinates of chlorine it would be 
necessary to choose the initial values as accurately as possible. The x co- 
ordinate was obtained from the position of the Cl—Cl peak in the c-axis 
Patterson projection. The x, z co-ordinates could have been similarly 
obtained from the b-axis Patterson projection, but, as mentioned earlier, 
the intensity data for this zone were not ready and therefore the correspond- 
ing Patterson map was not available. Now it was thought that instead of 
calculating the b-axis Patterson at a later stage to get the z co-ordinate of 
chlorine one could as well use the z co-ordinate of bromine as they would not 
be very much different from each other. Thus the x co-ordinate was found 
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to be 0-418 (from the c-axis Patterson) and z was assumed to be 0-437 which 
corresponded to the final refined value of bromine co-ordinates (Part I). 
The x and z co-ordinates of the light atom were taken from the hydrobromide 
at a stage when R(A0/) was 15-3 per cent. 


These set of initial co-ordinates (Table I) gave R® (AO/) = 37 per cent. 
for the structure factors F®) (h0/). The R value was rather large compared 
with that of 15-3 per cent. for the hydrobromide corresponding to the same 
light atom co-ordinates. The value was large even if we allowed for the 
fact that the hydrobromide was a compound containing a heavier atom. 
From this it could be inferred that either all the light atom positions were 
somewhat different from those assumed from the hydrobromide or that the 
co-ordinates of chlorine themselves might be appreciably different. The 
latter was thought more probable since the co-ordinate (z) of this atom was 
not obtained directly but was assumed to be same as that of bromine. In 
order to clarify the position a difference Fourier synthesis (Fig. 1) was com- 
puted using the structure factors F™) (A0/). As expected, it showed a large 
slope at the position of chlorine, the shift indicated being mainly in the 
z direction. Though slopes at the positions of the other light atoms were 
also noticed, it was thought that it would be better to refine the position of 
chlorine first, since it was quite heavy and then take up the refinement of the 
light atoms. 


The least squares formula, 


Fe 
me ZW? (Fo — Fe) | (5) a 


wey 
Ww =) 


was used to refine the x, z co-ordinates. In using the above formula the 
precautions mentioned in connection with the hydrobromide (Part I) were 
observed. 


Three cycles of least squares refinement of the position of chlorine alone 
were carried out. The R values at these stages were respectively R‘®) (h0/) 
= 33-6 per cent., R‘ (h0/) = 27:7 per cent. and R“) (h0/) = 27-2 per cent. 
The total shifts in the z and x co-ordinates were respectively 0-16A and 
0:06 A. Since the refinement of the co-ordinates of chlorine seemed to 
have converged that of the light atoms was taken up. 
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TABLE II 


of L-tyrosine hydrochloride 


Co-ordinates at different stages of refinement of the [001] projection 












Atom 1 2 3 4 5 

Cc, 0-030 0-035 0-035 0-035 0-035 

Cc, 0-130 0-127 0-127 0-127. 30-127 

Cs 0-130 0-128 0-128 0-128 0-128 

C, 0-030 0-032 0-032 0-035 0-035 

C, 0-930 0-935 0-936 0-936 0-941 

C, 0-930 0-940 0-936 0-936 0-940 

C, 0-030 0-035 0-035 0-035 0-035 Fractional 
C, 0-090 0-094 0-088 0-087 0-087 y co-ordinates, 
Cc 0-250 0-240 0-244 0-247 0-247 (y/b) 

O, 0-030 0-032 0-034 0-034 0-034 

O, 0-345 0-349 0-354 0°354 0-353 

QO, 0-275 0-278 0-282 0-278 0-273 

N 0-070 0-077 0-074 0-074 0-074 

Cl 0-250 

Cc. 0-098 0-100 0-100 0-100 0-100 

Cc 0-083 0-082 0-080 0-080 0-080 

C; 0-977 0-974 0-971 0-971 0-971 

Cc. 0-879 0-877 0-877 0-876 0-878 

c. 0-900 0-896 0-896 0-897 0-897 

& 0-003 0-008 0-010 0-010 0-010 

Cc. 0-763 0-765 0-767 0-763 0-763 Fractional 
C, 0-667 0-663 0-657 0-654 0-651 x co-ordinates, 


(x/a) 














ae ea” eo lu 








The Crystal Structure of L-Tyrosine Hydrohalides—II 25 


REFINEMENT OF THE CO-ORDINATES OF THE LIGHT ATOMS 


The first cycle of least squares refinement of the light atoms was carried 
out and it gave considerable shifts in all the light atom co-ordinates (Table I). 
A significant feature noticed in this connection was that the z co-ordinate 
shifts in all the light atoms were more or less of the same magnitude and 
also for most of them the sense was the same. This indicated that the mole- 
cule had to be shifted as a whole from the initial assumed position. The 
maximum shifts were about 0-16 A and 0-18A respectively in the z and 
x co-ordinates. These refined co-ordinates resulted in a reduction of the 
R value from R (h0/) = 27:2 per cent. to R® (hO/) = 24-5 per cent. 














Fic. 1. The b-axis difference Fourier projection using the structure factors F“ (h0/). Contours 
are at intervals of 0-5 e/A?. 


A second cycle of refinement of the light atoms was now gone through. 
The maximum shifts were found to be about 0-13 and 0-11 A respectively 
in the z and x co-ordinates. The structure factors F (h0/) were now cal- 
culated and the value of R‘ (40/) was found to be 22-8 per cent. 


In both the above stages of light atom refinement the overlapping group 
of atoms C,, Cy and O, was not included but was reserved for a later stage. 


Approximate positions of hydrogen atoms were now fixed up from bond 
length considerations. When these were included in the calculation of the 
structure factors, the R value obtained was R, (40/) = 21-9 per cent. 


Since the convergence seemed to have slowed down it was now decided 
to include the low angle reflexions also in the refinement. This would be 
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possible now since the hydrogen contributions had also been included. The 


overlapping group of atoms C,, Cy, O, was also now included in this third 
cycle of refinement. 


The method of treating the overlapping group was the same as that 
described for the hydrobromide compound (Part I). The sixth order normal 
equations of least squares involving the x and y parameter shifts of the three 
atoms C,, C, and O, were formed and were solved by the Gauss-Seidel itera- 
tive method (Whittaker and Robinson, 1945). The light atoms were also 
refined as usual and the new set of co-ordinates was used to calculate the 
structure factors F” (A0/). In spite of the fact that some of the shifts were 


quite considerable the value of R rather surprisingly remained stationary at 
R™ (hO/) = 23 per cent. 


It was now noticed that at high angles the calculated structure factors 
were systematically larger than the observed ones which indicated that the 
temperature factor needed alteration. Till now, only an isotropic tempera- 
ture factor with B = 2-35 A? had been used. Since the refinement seemed 
to have reached a fairly advanced stage the structure factors F‘” (A0/) 
were used along with the values of | F,| in plotting the usual graph 
of log {(| Fe |)/{| Fo |)} against (sin?0/A*). The graph gave a value of 
3-00 A? for B. The scale factor was also found to be altered slightly and 
when these were introduced into the new calculations the R values were 
found to be R'® (hO/) = 22-3 per cent. and Ry‘ (40/) = 20-5 per cent. 


It might be remembered that the refinement of chlorine had been stopped 
(after three cycles of least squares refinement) at a stage when R value was 
27-2 per cent. The fourth cycle of refinement of chlorine alone was now 
carried out using the structure factors F,,‘% (h0/). The shifts obtained 
were quite negligible, namely, 0-005 A and 0-001 A respectively in the z and 
x co-ordinates. The refinement of chlorine position was therefore assumed 


to be complete. As the shifts were negligible no new calculations were made 
with the new co-ordinates of chlorine. 


It was now thought that, since the atoms C,, Cy and O, were fairly weil 
resolved in the other c-axis projection it would be better to use, for these 
atoms alone, the values of x refined from that projection. The refinement 
of the c-axis projection had actually reached a well advanced stage after two 
cycles of least squares refinement, the R values then being R‘? (hk0) = 18-2 
per cent. and R,‘* (hkO) = 16-8 per cent. (see next section). The x co- 
ordinates of C,, Cy and O, were thus assumed to be those used in the cal- 
culation of F‘ (Ak0). All other co-ordinates were the same as those used 
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for F® (h0/). The new calculation gave R® (A0/) = 21-0 per cent. and 
Rx® (A0/) = 19-2 per cent. The improvement in the R value for the struc- 
ture factors F®) (h0/) as compared with F‘” (h0/) showed that the x co-ordi- 
nates of C,, Cy and O, were now definitely more accurate than those used 
for F™ (h0/). 

However, even now the R value seemed to be rather large (21 per cent.). 
But one striking fact was that the inclusion of hydrogen contribution 
improved the R value nearly by 2-0 per cent. This seemed to be due to the 
use of the structure factors Fy, (40/) for refinement. It also points to the 
importance of including hydrogen contributions at quite an early stage of 
refinement. The above situation may be compared with that in the hydro- 
bromide (Part I) when even at the final stage of refinement (R ~ 12 per cent.) 
the inclusion of hydrogen contribution practically left the R value unaffected. 


The rather large value of R (19-2 per cent.) was actually found to be 
due to the presence of extinction effect in some of the reflexions. The graphi- 
cal method of correction (Part I) was adopted and the value of G was found 
to be 0-15x10-*. The value of G was used to correct the observed structure 
amplitude and it resulted in Ry®® (A0/) = 16-0 per cent. 


A fourth cycle of least squares refinement of the light atoms including 
the overlapping group was now gone through. However, while solving the 
simultaneous equations of the overlapping group the shifts in the x co-ordi- 
nates of C,, Cy and O, were assumed to be zero and this was because the 
refined values from the other projection were used in the calculations. The 
new calculations gave practically a stationary value of Ry (40/) = 16-0 per 
cent. 


DIFFERENCE FOURIER SYNTHESES AND TEMPERATURE FACTOR REFINEMENT 


A difference Fourier synthesis (Fig. 2) was now computed using the 
structure factors F, (0/). The reason for using the structure factors 
which included the hydrogen contributions was that the synthesis was in- 
tended for a study of anisotropy of thermal vibrations of the atoms. The 
synthesis proved to be really useful. The main features of the map were 
the following: 


(1) There was a marked anisotropy in the vibration of the chlorine atom. 


(2) A similar anisotropy was noticed for the light atoms also, more 
especially for the oxygens. 


(3) The shapes of the positive and negative peaks indicated that the 
directions of vibration of the light atoms were more or less the same as that 
of chlorine. 
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(4) The slopes at the positions of some of the light atoms had not 
vanished indicating thereby that those co-ordinates needed shifting. 











Fic. 2. The b-axis difference Fourier projection using the structure factors F,,"® (/0/). Contours 
are at intervals of 0-2 e/A?. 


The third of the features mentioned above seems to be rather significant. 
The fact that the directions of vibration of the different atoms are almost 
the same indicates that the whole process can be approximated to one of 
vibration of the molecule as a whole. This goes to justify at least partially 
the view (Wells, 1956) that an overall anisotropic temperature factor could 
be applied for a particular reflexion. Moreover the assumption of aniso- 
tropic parameters for individual atoms and their refinement demands an 
abundance of data which could be had only if they are collected three dimen- 
sionally. Therefore only the overall anisotropic thermal parameters B,,, 
Bis, etc., were assumed so that we could write. 


F, (hkl) = Fe’ (hkl) exp-(B,; h? + By: k? + Byg + Bish!) (2) 


where F,’ (hk/) is the calculated structure factor without applying the tem- 
perature factor. The parameters B,,, By, etc., were now refined by the 
method of least squares. The initial starting values of B,,, B,s, etc., were 
obtained from the isotropic temperature factor with B = 3-0 A? (see Part I). 


However, just for the sake of interest the anisotropic thermal para- 
meters b,,, 533 and 5,3 for chlorine alone were refined from this b-axis projec- 
tion. Surprisingly the value of 5,,, b3; and 5,3 obtained for chlorine agreed 
closely with the overall thermal parameters B,,, B;; and B,; (Table III). This, 
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of course, may be explained as being due to the fact that chlorine, being 
rather heavy, contributes predominantly to the structure factors. However, 
the individual anisotropic parameters for chlorine were never used in the 
calculations but only the overall parameters were used. 


TABLE III 


L-tyrosine hydrochloride : Overall anisotropic thermal parameters 
obtained by the method of least squares 














(A)? 
Zone* 
By Bop Bs Bis 
[001] me a .. 00-0052 0-00954 
[010] i n .. 0-0050 oe 0-0344 — 0-00370 
Value from isotropic B = 3-0 A? 0-00617 0-00919 0-0300 — 0-00033 
Value for chlorine... .. 0-00320 ae 0-0332 — 0-00350 





* The anisotropic thermal parameters for chlorine alone were obtained by the method of 
least squares for the (010) zone. However these values were not used in calculations. 


Before applying the anisotropic temperature factor it was thought that 
the difference map itself could be used to obtain the shifts in the various co- 
ordinates. The formula 


a [? tro — 0] 


x 
2" Po 

a9) i 
was used for this purpose. The usual approximation of d%p9/)x* = 2pp, 
was made and the values of p for different atoms were obtained graphically 


using the electron density distribution near the peaks (Part I). The values 
were p = 2°65, 2-8, 3:0 and 4-3 respectively for C, N, O and Cl. 


The new co-ordinates thus obtained were used in the calculation of the 
structure factors F“® (h0/). However, to see how these new co-ordinates 
alone affected the R value the original isotropic temperature factor with 
B = 3-0 A? was applied to the structure factors F® (40/) and the values of 
the reliability indices were R® (A0/) = 17-7 per cent. and Ry” (0/) = 
15-8 per cent. When the overall anisotropic factor was applied the values 
were R (40/) = 16-4 per cent. and Ry?” (h0/) = 14-4 per cent. 


(3) 





Ax; = 
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The fact that the inclusion of hydrogen contributions improved the 
R value nearly by 2 per cent. indicated that the difference Fourier could now 
be employed to locate these atoms. Accordingly a difference Fourier using 
F‘°) (A0/) was computed (Fig. 3). As expected, it showed peaks at assumed 
positions of hydrogens. The map also showed that the anisotropy had been 
practically corrected for, since the characteristic positive and negative peaks 








Fic. 3. The b-axis difference Fourier projection using the structure factors F%°) (h0/). The 


assumed positions of hydrogen atoms (Table V) have been marked. Contours are at 
intervals of 0-2 e/A?. 


had almost completely vanished. The map gave slight shifts in some of the 
light atom co-ordinates. The new set of calculations gave almost stationary 
values of R™ (h0/) = 16-2 per cent. and Ry (40/) = 14-3 per cent. The 
refinement of this projection was now assumed to be complete. The values 
of R when reflexions with |Foi| = 0 are omitted are: R (h0/) = 14-5 per cent. 
and R,, (A0/) = 13-0 per cent. 


REFINEMENT OF THE c-AXIS PROJECTION 


In order to reduce the total work involved in this non-centrosymmetric 
projection it was decided to start the refinement with a fairly close approxi- 
mation to the real co-ordinates. This was possible since the refinement of 
the other projection had been started earlier and had reached a fairly well 
advanced stage. The x co-ordinates of the light atoms were now assumed 
to be those used in the calculation of the structure factors F’” (h0/). The 
y co-ordinates were chosen from bond length considerations. This set of 
initial x, y co-ordinates (Table II) was used to calculate the structure factors 
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F® (hk0). As mentioned in the last section the structure factors F‘ (h0/) 
had been used along with the values of |F,| to obtain more accurate scale 
and temperature factors. A similar graphical procedure was now adopted 
for this c-axis projection also. The value of B = 3-0 A? was obtained from 
the graph and this was indentical with the value obtained from the other 
projection. 


Using this isotropic temperature factor the R value was found to be 
R® (AkO) = 23-0 per cent. 


The least squares refinement of both x and y co-ordinates from this 
projection was now started. For the refinement of the y co-ordinates from 
this non-centrosymmetric projection the formula is the same as that given 
by equations (8) ¢nd (9) of Part I. However, for the refinement of the x co- 
ordinates, which was not carried out in the case of the hydrobromide, we 
have to use the more general expression 


DPR a ts +o) ‘ 
dD, Ge) are my 


dA ; 
(x) = 4th f, cos 2xky; sin 27hx; 





Ax; = 


where, for k even, 


)= = — 4nh f; sin 2xky; sin 27hx; 


and for k odd (7) 
dA . 
(3) = — 4rh fj; sin 2xky; cos 27hx; 


dB 

om) — Anh f; cos Inky; cos 2nhx; 

The refinement of y co-ordinates is thus easier since the quantities A, Aj, 
B, B; in (4) would already have been prepared during structure factor cal- 
culations while, for the refinement of the x co-ordinates, the quantities 
(dA/dx;), (2B/dx;) in (6) have to be calculated afresh using (7). 


The set of x, y co-ordinates obtained from the first cycle of least squares 
refinement was used to calculate the structure factors F‘®) (hkO) and it 
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resulted in R‘®) (AkO) = 19-0 per cent. When hydrogen contributions 
were now included, it gave R,,‘”) (AkO) = 18-4 per cent. A second cycle of 
least squares refinement now reduced the R value to F® (AkO) = 18-2 per 
cent. and R,* (Ak0) = 16-8 per cent. 


Some of the reflexions which showed extinction effects were now cor- 
rected for by the graphical method mentioned already. The value of G 
was found to be 0:2x10-%. After extinction correction the R values were 
R) (hkO) = 15-4 per cent. and Ry” (hAkO) = 13-6 per cent. 


DIFFERENCE FOURIER SYNTHESES AND TEMPERATURE FACTOR REFINEMENT 


A difference Fourier synthesis (Fig. 4) was now computed to study the 
general anisotropy of thermal vibrations. The map showed anisotropy for 
chlorine and also for some of the light atoms. This, of course, was expected 




















| | 


Fic. 4. The c-axis difference Fourier projection using the structure factors F® (hk0). 
Contours are at intervals of 0-2 e/A?. 


since a difference Fourier of the other projection at a corresponding stage 
also showed similar features. In this c-axis projection also the general nature 
of the vibrations of the light atoms and of the chlorine atoms were quite 
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similar. Co-ordinate shifts were also obtained by the usual method. How- 
ever, the use of the new co-ordinates left the R values practically unaltered 
at R® (440) = 15-3 per cent. and Ry) (kO) = 13-8 per cent. Only iso- 
tropic temperature factor with B = 3-0 A? was used in the above calculation. 


The refinement of the overall anisotropic factor was then carried out. 
For this projection, the appropriate form of the factor was exp-(B,, h? + 
B,.k*). The initial values were obtained as usual from the isotropic value, 
B=3-0A?. The value of B,, refined from this projection agreed well with 
that from the other projection (Table II). Application of the overall tem- 
perature factor gave R‘®) (hkO) = 14-4 per cent. and R,,“*% (AKO) = 13-1 per 
cent. 

It could be noticed that, in this projection also, the inclusion of hydrogen 
contribution improved the R value nearly by 1-3 per cent. Therefore, with 
the object of locating the hydrogen positions, a difference Fourier was com- 
puted (Fig. 5) using F (hk0). However, probably because this projection 
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Fic. 5. The c-axis difference Fourier projection using the structure factors F“® (Ak0). Contours 
are at intervals of 0-2 e/A?. 
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was non-centrosymmetric, the map was not as good as that of the b-axis 
projection. Yet, faint maxima were discernible at the positions of hydrogen 
atoms. The anisotropy was also found to have practically vanished. Shifts 
in the co-ordinates of some of the atoms were also obtained. The new cal- 
culation gave R“ (hkO) = 14-1 per cent. and Ry‘ (hkO) = 12-9 per cent. 
The almost stationary value of R was taken as an indication of the conver- 
gence of the refinement process which was therefore stopped. The co-ordi- 
nates at different stages of refinement of this projection are listed in Table IT. 


It was now noticed that two reflexions in this zone namely, 230 and 330, 
had rather large discrepancies between their calculated and observed structure 
amplitudes. A careful examination showed that both of them occurred at 
low angles and also for both of them |F,| was greater than |Fo|. This showed 
that the effects of extinction were present in them. The reason why they 
had been left out during the previous graphical correction was found to be 
the following. In order to look for the reflexions that suffer from extinction 
we have to take the intensities which are not corrected for L-P factors, and 
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Fic. 6. The final c-axis Fourier projection of L-tyrosine hydrochloride. Contours are at inter- 
vals of 1 e/A?, starting from 2-0e/A?. For chlorine the interval is 2-0 e/A2. 
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not the intensities as judged from the tables of structure factors, |F | and 
|F,|. While the values might not be large if taken from tables of |Fo| and 
|F,| they might be quite large when L-P factors are taken into account. This 
was actually found to be the case with the reflexions 230 and 330. These 
reflexions also fitted well into the extinction graph and the original value of G 
could now be retained and applied to these reflextions also. This gave the 
values of R®® (hkO) = 13-45 per cent. and Ry,“ (2k0) = 12-25 per cent. 
When reflexions with |F,| = 0 were omitted the values were: R (Ak0) = 
11-5 per cent. and Ry (Ak0) = 10-3 per cent. Thus the final R values in 
percentage are, 





With Without 
hydrogens hydrogens 





R (402) All reflexions 14-3 16-2 
|F,| = 0, omitted 13-0 14-5 

R (hk0) All reflexions 12-25 13-45 
|F,| = 0, omitted 10-30 11-50 





The final set of accepted co-ordinates is given in Table [V. The assumed 
co-ordinates of hydrogen atoms are listed in Table V. The calculated and 
observed structure factors are given in Tables VI and VII. The final Fourier 
syntheses are shown in Figs. 6 and 7. 


ACCURACY OF THE STRUCTURE DETERMINATION 


The standard deviations in the co-ordinates and in other parameters 
were calculated on the same lines as described in Part I. The values obtained 
are listed in Table VIII. The two values obtained for the standard deviation 
in the x co-ordinate which is common to both the projections agreed closely 
and only the average is given in Table VIII. It can be seen from the table 
that the standard errors in the different x, y, z co-ordinates for any particular 
atom are almost the same, so that an average value can be taken. Using 
these values, the standard errors in the various bond distances are found to be 


C—C = 0:033A, C—N=0-029A, C—O=0-027A 
and 
O—O = 0-021 A. 
The standard error in electron density is about 0-2 e/A? while that in a typi- 


cal angle in the benzene ring is about 2-3°. All the probable errors have 
to be taken as slightly larger for the overlapping group of atoms C,, C, and O,,. 


A4 
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TABLE [V 


Final accepted co-ordinates in L-tyrosine hydrochloride 








Atom x/a* y/b zie x (A) y (A) z (A) 
C, .. 0-0990 0-035 0-362 1-096 0-316 1-843 
C, .. 0-0815 0-127 0-157 0-902 1-147 0-799 
C; .. 09730 0-128 0-019 10-771 1-156 0-097 
C, .. 08785 0-033 0-088 9-725 0-316 0-448 
CG, .. 0-8970 0-941 0-304 9-930 8-497 1-547 
Cc, .. 0-0080 0-940 0-445 0-089 8-488 2-265 
C, .. 07620 0-035 0-922 8-435 0-316 4-693 
C, .. 0-6530 0-087 0-075 7-229 0-786 0-382 
Cc, .. 0-6690 0-247 0-155 7-406 2-230 0-789 
O, .. 0-2020 0-034 0-520 2-236 0-307 2-647 
O, .. 0+6390 0-353 0-003 7-074 3-188 0-015 
O, .. 07230 0-273 0-373 8-004 2-465 1-899 
N .. 0°5400 0-074 0-912 5-978 0-668 4-642 
Cl .. 0-4130 0-250 0-405 4-572 2-258 2-061 





* Average value from the two projections (Tables I and II). 


It is to be noted that the x co-ordinates refined from the two projections 
agree closely within the limits of the estimated errors. The final table of co- 
ordinates contains only the average value. 


It can be noticed that the final R values of the two projections are con- 
siderably different, namely, R,, (A0/) = 13-0 per cent. and Ry, (hk0) = 10-3 
per cent. The difference can be attributed to the fact that the R value in a 
non-centrosymmetric case is smaller than that in a centrosymmetric case. 
The statistical treatment of the problem (Luzzati, 1952) also indicates this 
result. Although Luzzati’s results are not strictly valid at the final stages of 
refinement and also for compounds containing unequal atoms, we can at 
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least take his results to be broadly true. The difference in the R values can 
be thought of as arising from the fact that, when the assumed structure is 
different from the correct one, only the differences in the modulii 
(\Fo| — |Fe|) can occur in a non-centrosymmetric structure while the vector 


TABLE V 


Assumed positions* of hydrogen atoms in L-tyrosine hydrochloride 








Atom x/a y/b z/c 

H, (O;) . 0°26 0-10 0-51 
H, (C,) 013 0-21 0-12 
Hg (Cs) re 0-96 0-21 0-90 
H, (C,) .. (77 0-13 0-88 
H;, (C,) .. 0°86 0-84 0-36 
H, (C,) 0-01 0-84 0-56 
H, (C,) 1. 0°74 0-97 0-80 
H, (C,) 063 0-01 0-13 
H, (Os) ae 0-75 0-36 0-46 
H,, (N) .. 0°50 0-15 0-05 
H,, (N) .. 0°50 0-15 0-73 
H,. (N) 0°56 0-99 0-82 





* (1) The positions of these hydrogen atoms are almost the same as in L-tyrosine hydro- 
bromide. 


(2) The position of Hy, given in this table was the one used in calculations and was originally 
assumed to form the bond N(M)H .... Clo9;._ However, it has since been assumed to form the 


bond N(M)H .... O,(Myo,;). No new calculations were made with this change in the position 
of this hydrogen. 


difference (F, — F,) occur in a centrosymmetric case. In the latter one it is 
obviously possible to have a finite |Fy — F,| while |F,| — |F,| = 0. Hence 
the observed mean value of | |Fo| — |Fe| | will be less for a non-centro- 
symmetric than for a centrosymmetric case. 
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TABLE VI—Contd. 














F F Fu F F Fu 
Observed Calculated Observed Calculated 
OS h06 
0 3-3 3°8 3°6 0 3-2 — 2:33 — 2:3 
l 3°3 — 4-6 — 4-7 ] 5+3 5°5 5°5 
2 2:1 — 1-2 — 1-2 2 1-4 — 1-1 — 1-4 
3 2:1 3-1 3°3 3 4-0 3-6 3°5 
4 2:1 — 356 — 3:6 4 0-0 — 1:7 — 1-7 
5 6°3 6:3 6°1 5 2-9 — 2:6 — 2:6 
6 2-7 — 44 — 4-4 
7 4-0 4-2 4-4 
8 1-4 1:6 = 1:5 
9 19 —22 — 22 





* 


These reflexions have been corrected for extinction effects. 


TABLE VII 


Observed and calculated structure amplitudes and calculated phase 
anglest for the hkO reflexions of L-tyrosine hydrochloride 














| Fol |F| | Fx| a? Fol |F| [Fx| [a 
Calculated Calculated 
Al0 h30 
l 17-6* 17-8 18-3 345 l 14-6 11-9 13-5 270 
2 14-5 12-7 15-5 149 2 31-4* 34-] 34-5 160 
3 32°7* 32-5 33-4 196 a 26:2* 29-0 29-1 352 
4 24-9 29-2 29-2 3 4 22:4* 24-6 25-4 199 
5 9-] 11-1 9-0 213 5 6:5 6°8 6:1 32 
6 14-7 13-6 13-7 99 6 8-6 9-7 9-5 36 
7 16-6 14:5 14-4 359 7 9-2 9-] 9-8 100 
8 11-3 8-4 9-0 170 8 11-2 10-3 10-0 18 
y 22-0 22-7 23-0 334 9 6:4 6-9 6°8 316 
10 7:6 7-9 8-0 ZiT 10 8-6 9-9 9-8 21 
ll 2:3 2:4 2:6 30 11 4-4 4-0 4-3 115 
12 5-4 5:5 5-5 22 12 2:8 2°8 2:8 158 
13 5:0 5:0 4-6 196 
14 3-2 3-4 3°4 11 





+t The phase angles correspond to F,. 
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42 


|Fo| 


h60 

0 18-9 
l 14-0 
2 10-0 
3 11-1 
4 8-1 
5 5:8 
6 6°8 
7 8-] 
8 4-0 
9 3°7 
10 a2 
11 3°5 
12 5-0 
h80 

0 8-2 
l 10-6 


Fl 


|Fu| 





|a°| 
Calculated 

20-9 19-3 169 
15-0 14-5 202 
10-7 10-0 231 
10-8 10°4 243 
6-4 7-0 292 
5:0 5:0 136 
6-1 6-5 7 
11-2 9-6 219 
4-0 4:0 18 
4-4 4-3 105 
3-0 2:9 196 
4-0 3-5 329 
6:0 6-0 204 
10-0 10:0 27 
10-0 9-9 138 
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|Fol |F| |Fx| 


|o°| 
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278 
313 
141 
343 
221 


177 
175 


162 

74 
148 
225 
346 
166 
339 





* These reflexions have been corrected for extinction effects. 


Fic. 7. The final b-axis]}Fourier projection of L-tyrosine hydrcchloride. 
vals of 1-0 e/A? starting from 1-0 e/A?. 








For chlorine the interval is 2-0 e/A?. 





Contours are at inter- 
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TABLE VIII 


L-tyrosine hydrochloride : Standard deviation in various parameters 

















(A) 
Atom — 
a (x) o (y) oa (z) Average 
Cc 0-022 0-023 0-023 0-023 
N 0-018 0-018 0-016 0-017 
O 0-014 0-016 0-015 0-015 
Cl .. 0-0055 “i 0-0065 0-0060 
Bond (/) a (1) 
C—C * 0-033 (A) 
C—N ‘ia 0-029 
C—O - 0-027 
O—O ‘e 0-021 
¢ in a typical angle in the benzene ring = 2°3°. o(p) = 0-2 e/A?. 


DISCUSSION OF THE STRUCTURE 


Intra-molecular features.—The interatomic distances and bond angles 
are listed in Table IX and are shown in Fig. 8. The inter-molecular distances 
including hydrogen bonds are given in Table X. Figure 9 shows the packing 
arrangement of the molecules in the unit cell. The notations used in the 
case of the hydrobromide (Part I) are used in the figures and discussions of 
the present case also. 


Table IX shows that the bond lengths in the benzene ring vary between 
1-350 A and 1-404A with an average of 1-:385A. This average value 
agrees closely with that (1-39 A) found in a recent accurate determination 
of the structure of benzoic acid (Sim, Robertson and Goodwin, 1955). How- 
ever, the range as well as the average of these bonds are smaller than those 
of the hydrobromide. The fact that the bonds are rather small on the average 
is most probably due to the following reason. It might have been noticed 
that during the refinement of the present compound the hydrogen contribu- 
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TABLE IX 


Inter-atomic distances and bond angles in L-tyrosine hydrochloride 





— 





Bond (/) Atoms Angle 
C,G 1-350 (A) C,—C,—C, 120° 6’ 
C,—C,; 1-372 C,—C,—-C, 121° 26’ 
C,—C, 1-395 C,—C,—C, 117° 42’ 
ca. 1-399 a a 120° 0’ 
C;—C, 1-404 C.—C,—C, 118° 12’ 
C.-C, | -399 C.—C,—C, 120° 42’ 
C,—C, 1-522 0,—C,—C, 123° 24’ 
C,-C, 1-530 0,—C,—-C, 116° 12’ 
C,—C, 1-510 C,—C,—C, 119° 12’ 
C,—O, 1-374 C;—C,—C, 122° 24’ 
C,—O, 1-270 C,—C,—C, 112° 45’ 
C,—O, 1-265 C,—C,—C, 110° 6’ 
C,—N 1-484 C,—C,—N 110° 48’ 

N—C,—C, 108° 27’ 
C,—C,—O, 122° 0’ 
C.--h 117° 36’ 
O,—C,—O, 120° 80’ 











tions were included in the calculation of structure factors at quite an early 
stage of refinement (when R value was about 20 per cent.). Since these struc- 
ture factors were used in the refinements they should lead to the correct co- 
ordinates which otherwise would have been effectively displaced towards the 
hydrogens, if the hydrogen contributions had not been included in the siruc- 
ture factors. In other words, the systematic refinement after including the 
hydrogen contributions seems to be the reason for an effective smaller average 
value of the bonds and conversely if cycles of refinement are not carried out 
after the inclusion of hydrogens (as was the case with the hydrobromide) 
the whole benzene ring would appear slightly larger. This fact seems to have 
been ignored in most structure determinations resulting in a rather larger 
yalue for the average bond in the benzene ring. 








Fi 
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Fic. 8. Observed values of (a) bond lengths (in A) and (6) bond angles in L-tyrosine hydro- 
chloride. 


The equation to the plane of the benzene ring was determined by the 
method of least squares, and is given by 


z+ 1-143 y — 0-570 x = 1-587. 


The largest deviation of 0-007 A from this plane occurs for the atom C, and 
the mean deviation (r.m.s.) is only 0-0045 A. The deviations of the various 
atoms from this plane are, 


C, = — 0-0049 C, = 0-0056 C; = 0-006 
C, = — 0-0068 C;=0-0006 and C,=0-0043 A 


This establishes the almost perfect-planarity of the benzene ring. The atoms 
O, and C, also lie almost in the same plane, their deviations being respectively 
+ 0-08 and — 0-075 A. The angles in the benzene ring are also quite normal 
with an average of 119° 40’. 


The bond C,—O, has a length of 1-374 A and is in good agreement with 
the values usually reported for a single bond between an aromatic carbon 
and an oxygen. Examples of such values are 1-36 A and 1-37A in B-resorci-_ 
nol (Robertson and Ubbelohde, 1938), 1-36A in oa-resorcinol (Robertson, 
1936), 1-38 A in glycyl-L-tyrosine hydrochloride (Smits and Wiebenga, 1953) 
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and 1-36 A in some additional complexes of quinol (Palin and Powell, 1947). 
In view of the general agreement in the features of the hydrobromide and the 
hydrochloride, it can be said safely that the rather large value of 1-42A 
found in the former compound (Part I) for the above bond is not real, since 
the corresponding value is only 1-374A in the present case. 


TABLE X 


Inter-molecular distances, including hydrogen bonds and bond angles 
in L-tyrosine hydrochloride 


Hydrogen bonds Hydrogen bond angles 





N(M’)H.... Cl’ »» 3°32(A) Cy—N (M’)—Clo; 137° 
N(M’)H .... Clon 331 C,—N (M’)—Cl’ 86° 
N (M’) H .... O. (Mig) 2-79 Cs—N (M’)—O, (Mio) 124° 
O, (M’) H .... O, (M) 2-50 C,—O, (M’)—O, (M) 111° 
O,(M)H.... Cl 3-11 C,—O, (M)—Cl 120° 
Short Contacts 

CI—Cy (Mig) 3-14 (A) 

N (M’) —Cloo, 3-39 

C1—Oz (Mio) 3-42 

CI—Og (Myo) 3-45 

Cy (M)—C; (M) 3-62 

C, (M)—O, (M) 3-51 





All the C—C distances are close to the standard value. The angles 


at the aliphatic carbons are also in agreement with the tetrahedral angle of 
109° 28’. 


The features of the carboxyl group are quite similar to those found in 
the hydrobromide (Part I). The Cy,—O, and C,—O, distances are 1-270 
and 1-265 A respectively and are in agreement with the values usually found. 
The slightly larger values for them as compared with the pure C=O distance 
(1-215, Pauling, 1945) are also compatible with the finding (see discussion 








later 
angl 
the | 





The Crystal Structure of L-Tyrosine Hydrohalides—II 47 






later) that both the oxygens are involved in hydrogen bond formation. The 
angle 121° for O.—C,—O, may seem to be rather small, at first sight. Now 
the standard error in bond angles is about 2°. Also, since two of the atoms 
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Fic 9, Unit cell of L-tyrosine hydrochloride projected down (a) the c-axis and (b) the b-axis. 
Hydrogen tonds are shown by dotted lines. 
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C, and O, are involved in overlap in one projection, the error in the angle 
can be taken to be slightly larger than 2° and the above deviation of about 
4° from the standard expected value of about 125° (Pauling, 1945) is not 
significant. Moreover, values of 122° in norleucine (Mathieson, 1953), 121° 
in DL-a-methionine (Mathieson, 1952), 122° in DL-8-methionine (Mathie- 
son, 1952) and 122° in glycyl-L-tryptophan dihydrate (Pasternak, 1956) have 
also been reported, which are all distinctly less than 125°. 


The equation to the plane of the carboxyl group is given by 
z + 0-020 y — 2-049 x = 8-269. 


The largest deviation from this plane is for the atom O, (0-032 A) and the 
r.m.s. deviation is 0-017A. The deviations of the various atoms are 
C, = —0-004 A, C,=0-032A, O,=—0-001A and O, = —0-017A. 
The nitrogen is out of this plane by 0-71 A. As the value seems to vary 
within wide limits in the various amino acids investigated so far, the value of 
0-71 A in our case does not seem to require any special attention. 


The aC—NH;>* distance is 1-484 A and is in agreement with that found 
in the hydrobromide (Part I) and other compounds. 


Molecular arrangement and hydrogen bonds.—As is to be expected, the 
molecular packing arrangement and the scheme of hydrogen bonds are, in 
the present case, identical with those found in the hydrobromide (Part I). 
They are therefore discussed only briefly here excepting where they require 
some explanation. 


The rather short OH .... O bond between the oxygen O,(M) of the 
phenolic hydroxyl group and the oxygen O;(M’) of the carboxyl group 
found in the hydrobromide (Part I) is remarkably confirmed in the present 
case. The value found now, namely. 2-50 A, is identical with that found in 
the hydrobromide. The O,(M)H .... (C1 bond is of length 3-11A and 
agrees well with that found for the hydrobromide namely 3-23 A, provided 
we allow for the difference in the ionic radii of bromine and chlorine, which 
is 0-1 A (Pauling, 1945). 


The angles C,—O, (M)—Cl (120°) and C,—O, (M’) — O, (M) (110°) are 
in good agreement with the tetrahedral angle. 


The hydrogen bonds involving the NH;* group are essentially the same 
as in the previous compound. The slight difficulty experienced in the case 
of the hydrobromide (Part I) in deciding between NH .... Oand NH .... Br 
bonds practically vanishes here. This is due to the fact that 
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N (M’)—O,(Myjpo,) distance is 2-79 A and considerably shorter than 3-01 A 
found in the hydrobromide. This only confirms our previous conclusion 
that the N(M’)H .... O,(Myp;) is the real hydrogen bond and N(M’)H 
...+ Clog, is not. In fact, the present value of 2:79 A for the NH .... O 
bond agrees closely with those observed in other compounds. Examples 
are: 2-79, 2-81 and 2-87 A in DL-serine (Shoemaker et al., 1950),2-73,2-87/and 
2:87 A in norleucine (Mathieson, 1953), 2-80 and 2:77 A in DL-methionine 
(Mathieson, 1952) and 2-75-2-89 A in glycyl L-typtophan dihydrate (Paster- 
nak, 1956). The slightly shorter length of NH .... O bond is also indicative 
of the bond being quite strong. The N—Cl distances are N (M’)—CI’ = 
3-32 A, N(M’)—Cloo:’ = 3:3 A and N(M’)—Cloo, = 3:39 A. The corres- 
ponding distances in the hydrobromide are 3-45, 3-45 and 3-50A and hence 
the two compounds closely agree in their relative distances also. The fact 
that N (M’)—Clpo; is slightly larger than tae other two which are almost equal 
(3-32 and 3-31 A) also supports our previous conclusion that the former is 
not a hydrogen bond. The hydrogen bond angles near the NH,* groups are 
also quite satisfactory. 


There is also a short Van der Waals contact between Cy (My) and Cl. 
Its value is 3-14 A. This short contact has also been observed in the hydro- 
bromide (Part I) where it has the value 3-30 A. The difference between the 
two is also just the difference in the ionic radii of bromine and chlorine. This 
short contact is thus established by way of confirmation, in the two compounds. 
Examples of short contacts have already been cited (Part [). In 
fact, the short OH .... O bond coupled with the above short C,—X (X 
stands for Cl or Br) contact seems to be an important feature in the two 
structures. It shows how, under conditions peculiar to a particular structure, 
some of the bonded and non-bonded distances can vary over fairly wide limits, 
while, of course, most of the other bonds are normal. 


The author wishes to thank Professor G. N. Ramachandran for valuable 
suggestions. He is grateful to the Government of India for the award of a 
Scholarship during the tenure of which the above work was carried out. 


REFERENCES 
Buerger, M. J. .. Acta Cryst., 1951, 4, 531. 
Kartha, G. and Ramachandran, G.N. IJbid., 1955, 8, 195. 
Luzzati, V. .. Ibid., 1952, 5, 802. 
Mathieson, A. McL. -. Wild, 1982. §. 332. 


ee .. Ibid., 1953, 6, 399. 
Palin, D. E. and Powell, H. M. .. J. Chem. Soc., 1947, 208. 








50 


Pasternak, R. A. 
Pauling, L. 


Robertson, J. M. 
———— and Ubbelohde, K. R. 


Shoemaker, D. P., Barieau, R. E., 
Donohue, J. and Lu, C. S. 


Sim, G. A., Robertson, J. M. and 
Goodwin, T. H. 


Smits, D. W. and Wiebenga, E. H. .. 


Srinivasan, R. 
Viervoll, H. and Ogrim, O. 
Wells, H. J. G. 


Whittaker, E. T. and Robinson, G. 


R. SRINIVASAN 


Acta Cryst., 1956, 9, 341. 


The Nature of the Chemical Bond, \thaca: Cornell Univ. 
Press, 1945. 


Proc. Roy. Soc., 1936, 157 A, 79. 
Ibid., 1938, 167A, 122. 
Acta Cryst., 1953, 6, 241. 


Ibid., 1955, 8, 157. 


Ibid., 1953, 4, 531. 

Proc. Ind. Acad. Sci., 1959, 49 A, 340. 

Acta Cryst., 1949, 2, 277. 

Ibid., 1956, 9, 709. 

The Calculus of Observations, 4th Edition (Blackie), 1945. 




















INFRA-RED SPECTRA OF FLUOROXYLENES 


By M. R. PADHYE* AND T. S. VARADARAJAN 
(Department of Chemical Technology, University of Bombay, Bombay-\9) 


Received April 15, 1959 
(Communicated by Dr. R. K. Asundi, F.A.Sc.) 


ABSTRACT 


The infra-red spectra of four fluoroxylenes (1, 2, 6-, 1, 2, 3-, 1, 3, 4- 
and 1, 2, 5-isomers) have been recorded in the 2-15 region on a Perkin- 
Elmer Model 21 spectrophotometer. From among the bands recorded, 
those which could be correlated to the fundamental vibration of the 
constituent benzene and methyl groups are discussed in detail and com- 
pared with particular reference to the type of substitutions. 


SEVERAL methyl benzenes were studied wit refehrence to their vibration 
spectra by Pitzer and Scott.1 Among them were the xylenes whose vibra- 
tion frequencies he assigned and correlated to benzene and methyl vibrations. 
Recently Wilmshurst and Bernstein? have studied toluene and m-xylene, 
in the light of fresh frequency assignments by Ingold et al.* and Mair and 
Hornig.* 


The infra-red spectra of several fluorinated benzenes have been studied 
by Nielsen et al.®»® The main approach of these workers has. been the 
classification of the frequencies into the different symmetry species. Several 
workers have studied the vibrations belonging to the methyl group in parti- 
cular. Among them may be mentioned Fox and Martin’ and Sheppard and 
Simpson.’ The latter workers classify the methyl vibrations as the internal 
and external vibrations. Randle and Whiffen? have studied a number of 
compounds with particular reference to the CH; rocking modes. Barcelo 
and Bellanato” in their study of methyl amines make a comparison of the 
relative intensities in the Raman and Infra-red spectra of the internal C-H 
stretching and bending vibrations of the methyl group. ars 


It is proposed to study in the present investigation, the infra-red spectra 
of four isomeric fluoroxylenes. They may be called 1, 2, 6-, 1, 2, 3-, 1, 3, 4 
and 1, 2, 5-fluoroxylenes and are represented diagrammatically. on next 
page. 


* Present address: Indian Institute of Technology, Bombay 18. 
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F F F F 
HC. () cus ( ous () Fa 
AF \/\cHs i CH, \ 

CH, 


1, 2, 6-Fluoroxylene 1, 2, 3-Fluoroxylene 1, 3, 4—Fluoroxylene 1, 2, 5-Fluoroxylene 


(I) (II) (III) (IV) 
EXPERIMENTAL 


The compounds were prepared in this laboratory from the correspond- 
ing xylidines and were found to be pure. The infra-red spectra were recorded 
on the Perkin-Elmer Model 21 double beam instrument with linear wave- 
length scale, principally using the sodium chloride prism. The C-H stretching 
frequency region was recorded also with the lithium fluoride prism. The 
spectra were recorded using the liquid as such with a compensating sodium 
chloride window in the other beam. Cell thicknesses of 0-1 mm. and 
0-025 mm. were used. In the case of very strong bands the substance was 
taken as a thin film in between two sodium chloride windows. The infra-red 
spectra in the sodium chloride region are reproduced in Fig. 1. 


RESULTS AND DISCUSSION 


Among the four compounds under consideration, compound I belongs 
to the C., group and the rest to C,. Out of the 48 fundamental vibrations 
for each molecule, 30 will belong to the phenyl part and 18 to the CH; pat. 
For the C,, group the 48 vibrations get divided into the respective species 
as 16a, +7a,+ 105,+ 105,. Here the x-axis is taken perpendicular to 
the plane of the molecule and the z-axis is taken as the two-fold axis as 
suggested by the Joint Commission on spectroscopy. All are allowed in 
the Raman effect and all but the a, vibrations in the infra-red. For the C, 
group the species division will be 31a’ + 17a’. 


The main interest in the present study of these spectra is to view them 
in the context of the analyses of the spectra of xylenes by Pitzer and Scott! 
and m-xylene by Wilmshurst and Bernstein.? The latter’s assignments are 
consistent with vapour band contours, depolarisation ratios, the Teller- 
Redlich product rule and a modified form of this rule suggested by Pitzer 
and Scott. Also of interest will be the effect of fluorine substitution on the 
phenyl and methyl vibration frequencies. Now it is proposed to analyse 
separately the phenyl and methyl vibrations. In Table I are presented the 
frequencies recorded for the four compounds along with their intensities. 
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Fic. 1. Infra-red Spectra of Fluoroxylenes. 
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TABLE I 


Infra-red frequencies of fluoroxylenes (cm.~') 








Compound I Compound II Compound III Compound IV 
662 m 667 w 696 m 701 vw 
721 s 703 s 720 s 721 s 
768 s 773 s 749 m 756 m 
817 s 776 s 765 s 806 s 
828 m 884 w 781s 861 s 
891 vw 898 s 860 s 928 s 
920 m 959 w 931 w 946 vw 
963 w 989 w 947 s 995 s 
987 w 1022 m 1004 s 1043 m 

1006 w 1062 s 1012 s 1116s 
1038 m 1089 s 1024 s 1144 m 
1081 s 1159 w 1105 s 1157 w 
1095 s 1195 s 1151 s 1205 w 
1157 m 1224 w 1188 s 1215 w 
1193 s 1241 s 1199 s 1259 s 
1220 s 1255 s 1253 s 1279 w 
1272 s 1258 s 1276 s 1297 w 
1312 w 1274 w 1391 m 1383 m 
1383 m 1285 w 1416 m 1420 s 
1445 s 1332 vw 1460 s 1447 s 
1475 s 1389 s 1499 s 1453 s 
1484 s 1458 s 1506 s 1462 s 
1600 m 1477 s 1603 s 1517 s 
1629 vw 1513 m 1631 s 1589 s 
1650 vw 1550 m 1658 w 1642 s 
1658 w 1570 s 1698 vw 1684 w 
1686 vw 1587 s 1731 w 1722 w 
1736 vw 1631 s 1821 vw 1725 w 
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TABLE I—Contd. 





Compound I Compound II Compound III Compound IV 








1777 vw 1655 vw 1865 w 1737 w 
1808 vw 1695 vw 1968 vw 1743 w 
1850 w 1712 vw 1995 vw 1808 w 
1923 w 1730 vw 2049 vw 1811 w 
1997 vw 1759 ww 2094 vw 1834 w 
2009 vw 1793 vw 2129 vw 1884 w 
2042 vw 1836 vw 2195 vw 1976 vw 
2091 vw 1877 vw 2221 vw 1997 vw 
2141 vw 1912 w 2249 vw 2011 vw 
2188 vw 1915 w 2395 vw 2040 vw 
2257 vw 1958 vw 2415 vw 2070 vw 
2296 vw 1997 vw 2438 vw 2100 vw 
2325 vw 2053 vw 2518 vw 2136 vw 
2352 vw 2138 vw 2556 vw 2185 vw 
2380 vw 2150 vw 2631 vw 2232 vw 
2409 vw 2192 vw 2659 vw 2262 vw 
2432 vw 2247 vw 2695 vw 2298 vw 
2456 vw 2282 vw 2735 w 2344 vw 
2518 vw 2311 vw 2762 vw 2363 vw 
2541 vw 2432 vw 2869 m 2409 vw 
2567 vw 2478 vw 2923 s 2435 vw 
2597 vw 2531 vw 2945 s 2468 vw 
2645 vw 2567 vw 2975 s 2525 vw 
2680 vw 2587 vw 3029 m 2567 vw 
2739 w 2634 vw 3058 w 2735 w 
2856 m 2735 w 3115 w 2865 m 
2923 s 2762 vw 2927 s 
2953 s 2869 m 2949 s 
2980 s 2927 s 2971 m 
3029 s 2949 s 3025 m 
3048 s 2975 m 3053 w 
3067 m 3034 m 3086 w 


3071 m 
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The phenyl vibrations—The proposed fundamental vibrations belong- 
ing to the phenyl part, together with their symmetry and the corresponding 
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benzene modes, are given in Table II. 


B CH 


y CH 


BCC 


y CC 


vy C—X 


TABLE II 


Fundamental vibrations of phenyl part (cm.-") 











Species Compounds 
cS €& Ba Benzene I II Ill IV 
a’ a, ag 992 768 773 781 806 
a’ b, bug 1310 ©4641312 «2613352 ~—=Oi«i jj 1297 
a’ a, = Vig, 1485 1484 1477 1506 1517 
a’ b, e,- Vig, 1485 ? ? ? 
a’ a, e,* Vea 1596 1629 1631 1631 1642 
a’ b, e,* Ven 1596 1600 1587 1603 1589 
a’ ay e,* V9q 1178 1157 1159 1151 = 1144 
a’ b, e,* V9p 1178 
a’ a, oa Vige 1033 
a’ be es" Vig, 1033 ia ais ” - 
a’ bs bu M45 1150 1095 1089 1105 1116 
a’ b, a8 vs 1326 ? ? ? ? 
a” b, au My 671 
a” as e, Viog 850 
a” b, e, Vip 850 3 - ws sie 
a” as e,* Vite 970 891 884 860 861 
a” b; e,* Via 970 920 898 928 
a” b, bg vs 985 963 959 9472? 946 
a’ a, bu Mp 1010 1006 989 1004 995 
a’ ay e,* Vea 606 662? 6672? 696? 701? 
a’ bs e,* Ven 606 - 7 us ‘ea 
a” b, bog v% 703 721 703 720 721 
a’ as e,* Vise 404 
a’ b, e,* Vigo 404 Se ie sa = 
C-CH, 1193 1195 1199 1205 
C-CH, 1272 1255 1276 1279 
C-F 1220 1241 1253 1259 
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Pitzer and Scott have assigned the benzene breathing vibration of 
992 cm.-! at about 1000cm.-! in toluene and m-xylene. Wilmshurst and 
Bernstein have assigned a value between 700-800 cm.-! in the case of 
m-xylene. These latter authors observe that the C-CH, stretching vibration 
and the ring vibration interact bringing down the ring vibration frequencies. 
In the case of fluoroxylenes the value of about 760 cm. for this vibration 
seems reasonable in view of a recent systematisation by Padhye and Vilad- 
kar. It is observed that in compound IV the strong band at 806cm.* 
is assigned to this vibration in preference to a weak one at 756cm.-' One 
component of the benzene ¥,, vibration at 1485 cm. should normally be 
about 1450-1500 cm.-! and the other about 1375-1400cm.', in substi- 
tuted benzenes. This region is crowded because of the presence of the CHs 
internal bending modes. The band in the region 1450-1500cm. gives 
a broad appearance in I and II whereas it is well separated in III and IV. 
There seems to be a superposition of three bands in this region. Two of 
these around 1450 cm. belong to the CH; asymmetric bending modes and 
the third one corresponds to one component of the 1485 cm. vibration of 
benzene. The compounds III and IV show additional strong bands at 
1416 cm.-! and 1420 cm.-! respectively. These could be an internal methyl 
bending mode or correspond to the other component of the benzene 
1485 cm. vibration. From intensity consideration they are tentatively 
assigned to the CH, internal deformation mode. 


In the 1600 cm. region there are two strong bands corresponding to 
the benzene 1596 cm.-! vibration. The separation of these bands in xylenes 
and disubstituted benzenes is low. For trisubstituted benzenes the separa- 
tion is shown to be rather high.12_ In this respect it is observed that II and 
IV on the one hand and the other two on the other show some similarity. 


The weak bands around 1300cm.-! have been assigned to correspond 
to the benzene v3, 1310 cm.-! vibration. The compound IV however does 
not show this band. 


Among the C—C bending vibrations, the v,., 1010cm.— vibration of 
benzene appears at around 1000cm.— in these compounds. It is weak or 
absent in the I and II compounds and strong in the other two. The other 
C-C bending modes that can be identified are discussed at the end. 


To make the assignment of the C-H stretching modes of the phenyl 
part more convenient, the spectra of xylenes in this region have been recorded 
with the lithium fluoride optics. They show four bands in this region 
whereas in fluoroxylenes, there are only three as shown in Table III and it 
seems that in each case a band near 3010 cm, is lost on fluorine substi- 


58 M. R. PADHYE AND T. S. VARADARAJAN 


tution. The four fluoroxylenes can be grouped in two types for the analysis 
of these vibrations; those having substitution at 1, 2, 3 and those at 1, 2, 4, 
Compound I and II fall into the first type and IT and IV compounds fall into 
the latter. The former have three bands near 3060, 3040 and 3030cm— 
whereas the latter have them near 3090, 3060 and 3040cm-! It thus 
seems that in addition to the C-H valence vibrations lost in xylenes, in the 
* 1, 2, 3” type of compounds, the other component of V5) 3080 cm.— frequency 
is changed to a C-X mode, whereas, in the ‘1, 2,4’ type of compounds 
one of v, 3047 cm. vibration is lost. 


TABLE LI 
Phenyl v C-H frequencies (cm.-') 











Compounds 
o-xylene m-xylene p-xylene 
[ II Ill IV 
3100 + 3096 3096 3067 3071 3115 3086 
3053 a 3044 3049 3048 du 3058 3053 


3016) 3025 3021 3029 3034 3029 3025 
30165 | Unresolved 3012 3003 7 fs 


The C-H inplane bending modes are expected to occur in the region 
1050-1300 cm.'_ In the case of fluoroxylenes, we expect only three vibra- 
tions. In all the compounds of the present investigation we get a band 
around 1100 cm.-! which may be assigned to the y,,, 1150cm. vibration 
of benzene. Again in all the compounds there is a band around 1150 cm. 
This can be correlated with vq, 1178cm.' vibration in benzene. It may 
be noted that the band at 1151 cm.‘ in compound III is quite strong when 
compared to the other three compounds. The third frequency belonging 
to this mode cannot be fixed because of the overlapping by the C-CH, and 
C-F stretching vibrations in the region. 


We expect three vibrations belonging to C-H out of plane mode. We 
get in all the cases one band between 860-890 cm.—! and another between 
900-930 cm.-' These may correspond to the two components of v,, 970 cm.-! 
vibration of benzene. There are weak bands at 963cm.-!, 959cm—! and 
946 cm.-* in compounds I, If and IV. These may correspond to the v,, 
985 cm. vibration of benzene. There is a strong band at 947 cm— in the 
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fl compound. In addition we have bands around 820cm.-? of medium 
intensity in compound I. 
TABLE IV 


CH, vibrations (cm.~') 











Mode Compound I Compound II Compound III Compound IV 
Asymmetric 2980 2975 2975 2971 
stretching 
2953 2949 2945 2949 
2923 2927 2923 2927 
Symmetric 2856 2869 2869 2865 
stretching 
Asymmetric 1475 1458 1460 1453 
bending 1445 ‘a 1416 1420 
Symmetric 1383 1389 1391 1383 
bending 
Rocking 1081 1062 1024 1043 
1038 1022 947 ? 
Torsion 





Among the C-X vibrations, where X is the substituent—in this case CH, 
and F—the bending ones are all below 600 cm.“ and so fall out of the range 
covered in the present investigation. The C-F stretching mode has been 
fixed at 1220, 1241, 1253 and 1259cm.—' in compounds I, II, DI and IV 
respectively. This is in accord with the values assigned by Nielsen et al.® ® 
for this vibration. 


Randle and Whiffen® attribute a frequency of 793 cm. in hexamethyl 
benzene to the C-CH, stretching mode. On the other hand Hastings and 
Nicholson!* have calculated this frequency to be 1385cm.-! In toluene, 
Pitzer and Scott' as well as Wilmshurst and Bernstein? assign a frequency 
around 1210 cm.~! to this mode. Again the latter workers assign 1170 cm.— 
and 1248 cm.—' to these modes in m-xylene. Consistent with these observa- 
tions, it is proposed in the present series of compounds to assign one fre- 
quency around 1190cm.—' and another around 1270cm.—! to these modes, 
The C-X stretching modes have also been included in Table I. 


60 M. R. PADHYE AND T. S. VARADARAJAN 


CH, internal vibrations—The assignment of fundamentals belonging to 
the CH; groups are shown in Table IV. There must be six stretching fre- 
quencies, two of symmetric and four of asymmetric type. Whereas Pitzer 
and Scott have kept all the frequencies coincident at 2950 cm.—!, Wilmshurst 
and Bernstein account for these with two vibrations at 2923 and 2953 cm-1 
In the present investigation four vibrations have been recorded in this 
region with the help of lithium fluoride optics, in all the cases. The vibra- 
tion at about 2860cm.-' may be taken as the symmetric stretching and the 
vibrations above 2900 cm.-' may be taken as the asymmetric stretching, 


The band at 2860 cm.-' corresponding to a polarised Raman line in 
m-xylene has been assigned as an overtone by Wilmshurst and Bernstein, 
and they have assigned the band at 2923 cm.~' to symmetric C-H_ stretching 
of the methyl group. The 2923 cm. band is polarised in the Raman spec- 
trum. Several earlier workers, on the other hand, have assigned a frequency 
around 2860 cm.—! to the symmetric C-H stretching of the CH, group and 
the bands above 2900cm.-' to asymmetric C-H stretching in a methyl 
substituted aromatic compound. 


The polarisation data in such cases is however tentative. Apart from 
the fact that there is a large amount of overlapping of the bands, the 
symmetry assignment of the internal modes of substituent groups have to 
be considered from the over-all symmetry point of view. This is specially 
important in non-planar substituents. In the present case one of the 
asymmetric C-H stretching vibrations of the methyl group under reduced 
symmetry would reduce to a totally symmetric vibration. Hence the vibra- 
tions around 2920 cm.-' have been assigned to the asymmetric C-H stretch- 
ing vibration of the CH; group. 


As for the bending modes we get in all cases a band around 1385 cm." 
which can be attributed to the symmetric bending mode. Again the bands 
around 1450cm.~ in all the compounds have been assigned to the unsym- 
metrical bending mode. We get two bands each in the I, II and IV com- 
pounds, corresponding to this mode. 


The bands at 1038 cm.—' and 1081 cm.“ in I are assigned to the rocking 
mode in line with similar assignment by Wilmshurst and Bernstein in 
m-xylene. The bands at 1022 and 1062cm.- in the II compound and at 
1,043 cm.-! in the IV compound have also been assigned to this mode. In 
the case of compound III there is a band at 1024cm.-! and another strong 
band at 947cm.-! It is to be considered whether both of these can be the 
rocking modes, as per Randle and Whiffen’s contension. According to 
these workers, for compounds with two adjacent methyl groups without 
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any other group adjacent to them, one should expect two bands at 995 cm? 
and 1022 cm.-! corresponding to the rocking modes. There are in the case 
of III, three bands at 947, 1004 and 1024cm.-! It is likely that the methyl 
rocking vibration and the vibration around 1000 cm.~ interact to give the 
anomalous intensity pattern for bands at 947 and 1004cm.~ in this com- 
pound. 


In the case of IV, the intensity distribution of the bands suggests a similar 
interaction. 


A few more observations can be made about the spectra of these com- 
pounds. They arise out of the consideration that the fluoroxylenes can be 
compared to the corresponding trimethyl benzenes taking into account the 
fact that the masses of fluorine and a methyl group are nearly the same as 
has been done by Ferguson ef a/."* in the case of fluorotoluene. The 1, 2, 3- 
trimethyl benzene has two strong bands in the region 650-900 cm. at 709 
and 766cm.-! The 1, 2, 4- trimethyl benzene, on the other hand, has three 
strong bands in the same region at 703, 806 and 873 cm. The compounds 
I and II which correspond to 1, 2, 3-trimethyl benzene have strong bands near 
660, 720 and 770cm.-! The compounds III and IV which correspond to 
1, 2, 4-trimethyl benzene have bands near 700, 720, 760 800 and 860 cm. 
The two bands near 760 and 800 cm. possibly show that the bands at 768 
and 773 cm.~! in compounds [ and II respectively may be the superposition 
of two bands which are only resolved in the ‘1, 2, 4” type of compounds. 
The higher frequency is assigned to the ring breathing vibration and the 
lower one probably is a y-C-H vibration. Compounds III and IV show an 
additional strong band at near 860 cm. characteristic of this type of sub- 
stitution. This is corroborated by the band at 878 cm.“ in 1, 2, 4-trimethyl 
benzene. 


Both types of trimethyl benzenes show only one strong band around 
700 cm.-! Among the xylenes only m-xylene shows a strong band in this 
region. In fluoroxylenes there are two bands in the region. One near 
660 cm.-! in 1, 2, 3 type of compounds or near 700 cm.“ in 1, 2, 4 type of 
compounds and the other near 720cm.™ in all the four compounds. It is 
likely that the band near 720 cm.—' may correspond to the band near 700 cm.-! 
in the 1, 2, 3- and 1, 2, 4-trimethyl benzenes and which from the comparison 
of xylenes arises due to meta-substitution. The 660-700 cm.-! band may be 
due to fluorine substitution or one of the components of the benzene 
606 cm.-! vibration which has been displaced to this region. 
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1. INTRODUCTION 


IN a nuclear magnetic resonance experiment, when the coil containing the 
sample is resonated by a capacitor and the circuit fed from a constant current 
source, the voltage drop across the coil is given by the approximate relation 


v = vo [1 — 140QfX], 


assuming |4rQfX| <1. Here vo is the voltage far off nuclear resonance, 
Q is the quality factor of the circuit, X is the complex r.f. nuclear suscepti- 
bility and f is the filling factor representing the fraction of the coil’s magnetic 
energy stored in the sample. We note that both the absorption and dis- 
persion components of the nuclear magnetic resonance signal are proportional 
to the Q of the circuit. 


Due to skin-effect, the dielectric losses in the coil form and the presence 
of a metallic shield round the coil, the values of Q attained by experimenters 
in this field have rarely exceeded 250, a common value being about 100. 


While improvements in the mechanical and electrical design of the coil 
are not likely to be of much help, it is possible to increase considerably the 
effective Q of the circuit by the use of an active network to supply energy 
to the tuned circuit to compensate for its losses. This can be achieved by 
introducing a shunt negative resistance in parallel with the coil. 


2. PRINCIPLE OF Q-MULTIPLICATION 


If we assume that the shunt resistance of the tuned circuit at resonance 
is R and that a negative resistance —R, is introduced in parallel with it, the 
effective shunt resistance becomes 


RR,» 
R, oe R ’ 





_ 


* Based on a paper read at the Low Energy Nuclear Physics Symposium, held at Bombay, 
on 27-28th February, 1958, under the auspices of the Department of Atomic Energy. 
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This corresponds to a multiplication of the original Q of the circuit by the 
ratio 


The Q-multiplication can in principle be made arbitrarily large by letting 
Rn approach R. If positive feedback is introduced through the resistor R; 
from an amplifier of gain A and internal resistance R; (see Fig. 1), it can be 
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Fic. 1. Basic circuit of the Q-multiplier using shunt negative resistance. 


easily shown that 


~~ Ry R; 
r= — A—1 A-—1?’ 
assuming the input impedance of the amplifier to be large compared to the 
other circuit resistances. As the circuit will break into oscillations when Ry, 
becomes equal to R, the gain of the amplifier has to be stabilised by a high 


degree of negative feedback (Ginzton, 1945) and a high —— resistor 
has to be used for Ry. 


3. NUCLEAR MAGNETIC RESONANCE EXPERIMENTS WITH A Q-MULTIPLIER 


Experiments on proton resonance in a sample of water, treated with 
ferric nitrate solution, were conducted with the amplitude bridge (Thomas 
and Huntoon, 1949). A Q-multiplier of the type described by Harris (1951), 
which employs a single tube acting as a cathode follower in the feedback 
patn, and consequently has a high inherent stability, was used. A minor 
modification to the Q-multiplier involved the application of feedback to the 
junction between two fixed capacitors connected in series across the main 
tuning condenser instead of splitting the latter. The complete circuit is 
given in Fig. 2. A brief analysis of the circuit, neglecting stray capacitances, 
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Fic. 2. Amplitude bridge with a Q-multiplier. 


shows that the effective negative shunt resistance introduced across the coil is 


_p = — 4. Reto t Re (% + 1) Re + tp} 
. +-h-, 


where rp and yu are respectively the plate resistance and amplification factor 
of the tube. 


The variable resistance R, served to maintain a given voltage across 
the coil for different effective values of Q. At a frequency of 7-35 Mc./S., 
a maximum increase in signal strength by a factor of about 5 was observed. 
The Q-multiplication achieved was limited because the shunt resistance of 
the tuned circuit (about 3,500 Ohms) was smaller than the magnitude of the 
minimum shunt negative resistance that could be generated by this simple 
Q-multiplier. It is) however, ideally suited for work in low field NMR 
spectroscopy, for the shunt resistance of the tuned circuit encountered there 
is quite high, a typical value being 780 K2 at 50 Kce./s. (Brown, 1950). At 
frequencies of the order of a few megacycles, it is desirable to use an amplifier 
with a stabilised gain of about 10. 


Another design criterion that has to be borne in mind at high frequencies 
is to satisfy the condition that the circuit be fed from a constant current 
source even at high values of the Q-multiplication factor for maximum benefit. 
A 0:1 mmf. capacity between the r.f. input point and the terminal of the 
tuned circuit at 20 Mc./s. corresponds to a reactance of about 80 K2 whereas 
the effective shunt resistance of the tuned circuit at a Qes; of 5,000 may be 
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of the order of 200 K2._ This difficulty does not arise in the nuclear induction 
method of resonance detection (Suryan, 1958). 


4. DISCUSSION 


Q-multiplication can enable effective Q’s of the order of 10,000 or more 
to be achieved without undue difficulty. This will lead to a hundred-fold 
or even greater increase in the strength of nuclear magnetic resonance signals. 


If q is the Q-multiplication ratio, it may appear from the formula of 
Bloembergen ef al. (1948) that, for a detecting system whose band width is 
determined by the final phase-sensitive detector, the signal-to-noise voltage 
ratio should also increase by a factor 1/q. This is, however, not so because 
both the nuclear magnetic resonance voltage and the thermal noise voltage 
generated in the sample coil are affected equally by the feedback. But the 
amplifier noise at the output will be somewhat less due to the Jower amplifi- 
cation now required. There might also be a slight improvement for an 
oscillographic display due to the increased selectivity of the circuit. 


Recently Benoit, Grivet and Guibé (1958) have independently appre- 
ciated the advantage of increasing Q by positive feedback in the design of 
a maser based purely on nuclear magnetic resonance. 


5. SUMMARY 


A method of increasing the strength of nuclear magnetic resonance 
signals by factors up to 100 or even more is suggested. The principle is to 
increase the quality factor Q of the coil containing the sample by means of 
positive feedback. Experiments on proton resonance in water, treated with 
ferric nitrate solution, using a Q-multiplier with the amplitude bridge, are 
described. 
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ABSTRACT 


Starting from the studies of Kleene and Mealy on sequential machines, 
in this paper is presented a formalism which, in a sense, unifies their treat- 
ments. From the specification of the required machine behaviour in terms 
of events and associated output states, a uniform procedure is given for 
obtaining a transition table and from that a minimal machine, whenever 
such a complete reduction is possible. The various steps of the synthesis 
procedure are so stated that they can be easily programmed on a computer. 


1. INTRODUCTION 


IN recent years, sequential machines have been extensively studied under the 
guise of nerve nets,’ switching networks? *;}* and automata.5 For the pre- 
sent purpose, we take as our point of departure the presentations of Kleene! 
and Mealy* and make free use of the concepts introduced by them wherever 
convenient. In several respects these two papers complement each other 
and, together, go a long way towards the solution of the synthesis problem, 
as we shall presently explain. 


A sequential machine is described in terms of a finite set of input states, 
a finite set of output states and a finite set of internal states. It is assumed 
that the present output state (O¢) and the next internal state (St,,) are uniquely 
determined by the present internal state (S;) and the present input state (I). 
Here, ¢ is the time parameter and it is understood that the machine is syn- 
chronous, in that, changes in the states take place at discrete instants of 
time #; (f= 1, 2, 3, ....). 


* Since we shall be concerned only with finite sequential machines in this paper, we shall 
henceforth drop the attribute ‘finite ’. 


When the first draft of this was ready, we received a copy of the paper “ Analysis of 
Sequential Machines—I ” by D.D. Aufenkamp and F. E. Hohn. Aufenkamp and Hobn are 
also concerned with the development of a formalism in terms of which the analysis of a sequential 
machine can be readily mechanised. However, their approach is entirely different from ours. We 
are grateful to Professor Hohn for making a copy of the paper available to us. 
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The problem of analysis of a given sequential machine is one of a 
complete description of its behaviour. This can be done by either of two 
equivalent methods: 


(1) The first method makes use of a Transition table or a state diagram 
which specifies exhaustively for each pair of present internal—present input 
states, what the associated pair of next internal—present output states is. 
This is the approach taken by Mealy. 


(2) The second method makes use of the notion of events. Events are 
input state sequences of finite length. For definiteness, assume that the 
sequential machine starts with a distinguished internal state. Then its com- 
plete behaviour is described by giving its output for each possible event. 
This is the method followed by Kleene. In his paper, he has answered in 
principle the basic question, what classes of events are representable by 
sequential machines? The formalism, developed with this in view, allows 
him to suggest a solution to the synthesis problem, viz., how to construct a 
sequential machine to represent a given event? However, his solution to 
this is mostly of academic interest. 


Neither of these two methods by itself seems to be adequate enough to 
answer satisfactorily the synthesis problem. In the synthesis of sequential 
machines one has three major aspects to contend with, viz., 


(a) specification of the required behaviour, 


(b) translation of the specification to a state diagram and the reduction 
of equivalent internal states, and 


(c) translation of the reduced state diagram to hardware. 


Mealy, in extending the works of Huffman and Moore, confines his 
consideration almost exclusively to aspects (5) and (c), i.e., given a transition 
table or, equivalently, a state diagram, he develops procedures to reduce it 
to its minimal form and to realise it in terms of circuitry. Kleene, on the 
other hand, is concerned primarily with aspect (a) of the problem, i.e., the 
development of adequate concepts to describe (and hence, equivalently, to 
specify) the machine behaviour in terms of its input states and output states, 


Thus, it would seem natural to look for a unified treatment which would 
exploit equally methods | and 2 described earlier and, hence, treat adequately 
all the 3 aspects of the synthesis problem. Clearly, a prerequisite for such 
a unification is the availability of a simple procedure which will enable one 
to go from a description of machine behaviour by method 2 to a description 
of its behaviour by method |. 
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In this paper we present a formalism which enables one to proceed from 
a specification of the required machine behaviour in terms of events and 
associated outputs to a transition table. From this initial transition table 
a reduction table is formed which, together with a reduction algorithm, allows 
one to arrive at the minimal equivalent machine whenever such a complete 
reduction is uniquely possible. We do not consider in this paper aspect (c) 


of the synthesis problem which calls for the realization of the reduced machine 
in terms of hardware. 


The formalism is developed such that the various steps in the synthesis 
procedure can be readily mechanised or programmed on a computer. This 
is an essential aspect of the problem in practice, since, without some such 
mechanical assistance, the synthesis of large, complex sequential machines 
will be almost impossible. 


2. CONVENTIONS 


Following the customary practice, we shall assume that the input, output 
and internal state organs of the machine are binary state devices, and refer 
to them as flip-flops. Thus, let the sequential machine have n input flip- 
flops and m output flip-flops. Then there are 2” possible input states and 
2™ possible output states. Call the former set I and index its members 
(Ip, I, .-.-, 14); call the latter set ® and index its members (4), 4), ¢., 


vee cy PyM_s). 


A k-sequence is a sequence of input states of length k. For a given k, 
every k sequence can be written as a k-tuple whose components read from 
left to right constitute the input states in that k-sequence. The 1-sequences 
are the individual input states. 


For each fixed k, the k-tuples are ordered according to the lexicographic 
ordering of their components. Thus, for each k, the first k-tuple is 


(Ih, I, .---, I) and the last (I,7_,, I,n_,, ...., 1). There are 2” 1-tuples, 
2” 2-tuples and in general, 2*" k-tuples. 
The integer N is a given parameter for the problem. N is the maximum 


length of an input sequence identified by the sequential machine as an event. 
Hence N is a measure of the finiteness of the memory of the machine. 


The set of input sequences of lengths k, for k = 1, .... N, has 


mn (>nN __ 
TWN) = HIM +... + IIHR (ALT D) 


members. Call this set Z and index its members with natural numbers 
(starting from 1). In each set of k-tuples, the indexing follows the ordering 
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defined earlier. The k-tuples themselves are ordered according to their 
lengths, starting with the 1-tuples and ending with the N-tuples. 


3. THe SPECIFICATION PROBLEM 


(a) Aspect.—It is assumed that the sequential machine always starts 
from a distinguished initial internal state denoted by S . The required 
behaviour of the machine is now specified by the assignment of 
outputs from the set ® to the elements i of J. Thus, if ¢q is assigned to i, it 
means that starting from the initial internal state S), the input sequence i 
terminates with the machine in the output state dg. 


In addition to the above assignments, certain end conditions will have 
to be prescribed. This is necessary because it was assumed that the 
sequential machine could identify k-sequences only fork <N. Hence, for 
completeness, its behaviour for k-sequences with k > N should be specified. 
(As we shall see, in practice, it is sufficient to specify for kK = N + 1 in addi- 
tion to K<N.) This is done, most naturally, by giving explicit rules which 
identify k-sequences, for k > N, with the elements i of 7 in a unique manner, 


Several cases will have to be distinguished according as the assignment 
of outputs is complete or partial. 


Case 1.—An output is assigned to every element i of J. This is the case 
of the ‘ exhaustive’ specification and admits of a complete reduction to a 
minimal machine uniquely. 


Case 2.—Outputs are assigned to some elements of J and it is specified 
that the assignment of outputs to the remaining elements of J are not of 
consequence (and, hence, can be done arbitrarily). The unassigned sequences 
are the so-called ‘“‘ Don’t-care ” sequences, and, clearly, any particular solu- 
tion will depend upon the way outputs are assigned to the ‘ Don’t-care’ 
sequences. Each such completion of assignment will reduce case 2 to an 
instance of case 1. 


Case 3.—Outputs are assigned to some elements of J and, in addition, 
auxiliary constraints are prescribed as holding between the elements of J, 
These constraints may be of the following type: 


(a) The elements i and k of J are indistinguishable or equivalent. By 
this is meant the following: The same output is assigned to i and k and, 
in addition, the subsequent behaviour of the machine is identical in the 
2 cases. Under these circumstances, it would be natural to say that i and k 
refer to the same event, 
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(5) Certain input sequences are forbidden, i.e., either they do not occu, 
as inputs to the machine or, if they do, the machine stops, say. 


(c) Other specialised constraints, if any. 


Clearly cases 1, 2 and 3 are not independent. In particular, the 
exhaustive specification of case | contains implicitly all the information 
required to partition the elements of J into equivalence classes as described 
in case 3(a). In fact, what we shall be concerned with in the rest of the 
paper is the determination of a uniform procedure to accomplish just this. 
Also it must be pointed out that where indistinguishability relations [as in 
3 (a)] are known to exist, the end conditions cannot be prescribed arbitrarily, 
but must be obviously consistent (see § 7 below). 


4. THE ASSIGNMENT OF INTERNAL STATES 


The specification of the required machine behaviour having been given, 
the crux of the remaining problem in synthesis is the assignment of internal 
states with their transitions such that the resultant machine is minimal. 
Minimality implies that no machine with a smaller number of internal states 
will be able to satisfy al] the specified requirements of behaviour. This 
problem is best handled in 2 distinct steps, as follows: 


Step 1.—First, a uniform procedure is given for assigning internal states 
to input sequences, i.e., to elements i of Z. 


Step 2.—Next, a transition table is formed and a uniform procedure is 
given for finding the equivalence among the internal states. This results in 
a table with a minimum number of internal states satisfying the specification. 


To begin with, we shall restrict our attention to case 1, i.e., that of the 
‘exhaustive’ specification, and develop the details of steps 1 and 2 for this 
case. Later, we shall show how these details can be suitably modified so as 
to reduce cases 2 and 3 to case l. 


5. CASE 1: ‘* EXHAUSTIVE’ SPECIFICATION 


Step 1: Assignment of internal states.—It is assumed here that an output 
from the set ® has been specified to each element i of Z. Now, to each ia 
distinct internal state S;,, is assigned, starting with S, for the l-tuple i = 0. 
S, is the distinguished initial internal state. Thus, Table I, the assignment 
table, is formed as shown in Fig. 1. 


Since the set has 


2n (m8 — | 
J(N) = — age ) 


members, Table I will have that many rows. 








the 
ion 
ded 
the 
his. 


ily, 


en, 
nal 
al. 
ites 


‘his 
ites 


> is 
; in 
on. 
the 
this 
) as 


put 
ia 


ent 








On the Synthesis of Finite Sequential Machines 





No. Input sequences Output states Internal states 
(elements i of 7) ; 








1 i=0 Pa, Si 
2 i=l be, S, 
3 i=2 $., Ss 
J (N) i=7J(N) — | baJin)-1 SZJin) 
Fic. 1, Table I: The Anslgnseumt Table. ‘ eniies 


Step 2: Construction of the transition table—Using Table I, we wish 
now to construct Table II, the transition table. The transition table, as 
shown in Fig. 2, is in the form of a matrix with 2" columns and J (N) + 1 
rows. The columns are named Ij, 1, ...., I,n_, and correspond to the input 
states with the same names. Similarly, the rows are named Sp, S,, ...., SZ 
and correspond to the internal states with the same names. Let the cell 
located at the intersection of row S; and column I; be named (Sj, I;). Each 
cell (S;, 1) contains two entries—one on the left, called the S-entry, which is 
the name of an internal state, and one on the right, called a ¢-entry, which 
is the name of an output state. 


The physical significance of the matrix is as follows: The names of the 
rows refer to the present internal states and those of the columns to the pre- 
sent input states. The S- and ¢-entries in the cell (S;, Ij) give the next internal 
state—present output state pair associated with the present internal state S; 
and the present input state Ij. 


Each row of Table I contains the names of one internal state and one 
output state. Each such pair S-¢ is entered as an S-¢ entry in one cell of 
Table II. To begin with, in the first row (named S,), the S-¢ pair corres- 
ponding to the 1-tuple j (in Table I) is entered in the cell (Sp, Ij). This is 
done for j= 0, 1, ...., 2*— 1. Any other cell (Sx, Ij) with k > 1 and 
j=0, 1, ...., 2"—1 is now filled in as follows: In Table I consider the 
input sequence i(k) which has been assigned the internal state S. Form 
the extended sequence by adding I; to the extreme right of i (kK). Let this 
input sequence be denoted by i(k) Ij. Now, fill in the cell (Sx, I), the S-¢ 
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pair in Table I corresponding to the sequence i(k)I;. In case i(k) is itself 
an N-tuple, i(k) I; will be a (N + 1)-tuple and, hence, there will be no row 
corresponding to it in Table I. However, by applying the proper end condi- 
tion to i (k) Ij, it will be identified with one of the rows of Table I uniquely. 
The S-¢ pair of that row is now used to fill in the cell (Sx, Ij). 





I,n_, 





Siva $; Sti Si p; 














Fic. 2. Table II: The Transition Table. 


The form of Table II shows that it is a transition table. The proof that 
the transition table so obtained does in fact satisfy the specification given 
follows directly from the explicit rule of construction of the input sequences 
which correspond to the internal states which are row names of Table II. 


Step 2(b): Construction of the reduction table-——Referring to Table II, the 
internal states which are the names of rows, are assigned superscripts accord- 
ing to the following rules. For definiteness, we shall use the Greek letters 
(in their alphabetical order) for this purpose. 


S; and S; have the same superscript if, and only if, the ¢-entries in 
those two rows, along corresponding columns, are identical. 


A uniform procedure for assigning the superscripts may be given. 
First partition the set of internal states (which are row names) according to 
their ¢-entries in the Oth column and assign superscripts. Next, re-partition 
these subsets according to their ¢-entries in the Ist column and assign fresh 
superscripts to the subsets newly formed and so on till the last column is 
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reached. The internal states which belong to the same subset, after the com- 
pletion of this operation, will have been assigned the same superscript. 


In the modified Table II the row names are internal states with super- 
scripts. From this modified Table II, the reduction table, Table III, is con- 
structed, as shown in Fig. 3. The reduction table is in the form of a matrix, 
with the rows corresponding to the next internal states (written down in order, 
starting with S,, together with their superscripts) and the columns corres- 
ponding to the input states written in order from left to right. In the cells 
are written down those internal states (with their superscripts) which associated 
with the input state which is their column name give rise to the next internal 
state which is their row name. 

. Qo a 
Nh | 
Sit ~*~ | 
\ 
| 





Io I, eee % In, 


S, Ss 








S,* 





S,7 








S?gin) 





Fic. 3. Table III: The Reduction Table. 


This can be done with the help of the modified Table II as follows: 
write down in cell (S;, 1;) of Table III all the row names of Table II whose 
S-entry in the column named I; is Sj. 


In the reduction process, we wish to merge all internal states which are 
equivalent to one another, into a single state. Thus when the reduction 
process is complete, we shall be left with a machine with a minimum number 
of internal states satisfying the specification. 


Clearly. the states (Sq,, Sa,, ...., Sa,) are equivalent if, and only if, they 
satisfy the 2 conditions: 


(1) For every fixed input state Ij, the pair (Sg, — Ij) results in the same 
Output state for all aj, i= 1, 2, ...., m. 


(2) For every fixed input state I;, the pairs (Sg, — Ij) result in the next 
internal states Sp, such that, either all the Sp, are the same, or are equivalent. 
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From the way the superscripts have been assigned it is evident that, to 
begin with, in Table III, all internal states satisfying condition 1 will have 
the same superscripts. However, in general, such states with the same super- 
scripts will not satisfy condition 2. The iteration procedure outlined below 


will enable us to modify the superscripts systematically, thus, completing the 
reduction process. 


In Table III, let the columns be ordered from left to right and in each 
column, the cells from top to bottom. We shall say a row has superscript a 
if the rowname is a next internal state with superscript a. 


The modification in the superscripts is done by the addition to them of 
indices i which are natural numbers. Thus, a re-partitioning of the set with 
superscript a; is done by introducing the superscript «;,,. Similarly for 8, 
y, etc. The superscripts are ordered according to the following scheme: 


’ B, Bi, Bs 


As each new superscript is introduced, it is added at once to this list of super- 
scripts at its proper place. 


The following iteration scheme is now used: 


(A) Find the first cell in column 0 containing an internal state with super- 
script a. Let its row have the superscript w. Find the next cell in the same 
column with an internal state also with the superscript a but whose row has 
a superscript different from w, say w,. Change all the superscripts 
a in that cell, and in all the succeeding cells, in the same column, whose rows 
have superscripts w,, to aj,, (where a up to a; are already in the list of super- 
scripts). Make corresponding changes in the superscripts of these internal 
states wherever they occur in Table’ III (including the names of rows). If 
no such next cell is found, of course, no modification is to be done. Repeat, 
now, the process with column |, then with column 2 and so on up to column 
2" — 1. 

Let A be called a minor cycle. 


(B) The minor cycle is repeated using a, in place of a and, when this 
is done, with the next superscript in the list in place of the previous one and 
so on till the entire list of superscripts is exhausted. Let A and B together 
be called a major cycle. 


(C) The major cycle is iterated till no modification of superscripts takes 
place in one complete major cycle. 


The reduction process is now complete and the internal states having 
the same superscripts are equivalent and hence can be merged. 
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The proof that this iteration procedure terminates and does result in 
complete reduction of the equivalent states is straightforward and we shall 
not go into the details here. In fact, the procedure given above is essentially 
a systematization of Mealy’s Rule I [pp. 1056 of (4)]. 


6. Case 2: SPECIFICATION WITH ‘ DON’T-CARE’ SEQUENCES 


Here, outputs are specified for some elements i of J and the remaining 
elements of J are specified as ‘ dont’-care’ input sequences. Table I is now 
constructed, as before, assigning a distinct internal state S;,, to each input 
sequence i. The ¢-entries for the ‘ don’t-care ’ sequences are written in as D. 


Thus, when Table II is constructed from Table I, some of the cells will 
contain D as their ¢-entries. Before forming the modified Table II, each 
entry marked D will have to be replaced by a particular 4%. Superscripts 
are now assigned to the rownames and the reduction process completed as 
before. Clearly, every such replacement of the D’s by the ¢’s will result in 
a reduced machine satisfying the given specification. 


The problem of great interest, in practice, is how to assign the ¢,’s to 
the D’s so that the reduced machine has a minimum number of internal 
states. Clearly, not much can be said about this in detail, since the solution 


for any specific case will depend on the way the D’s are distributed among 
the cells of Table II. However, in general, the aim should be to assign the 
¢,s So as to minimise the number of distinct superscripts in the modified 
Table I. 


7. Case 3: SPECIFICATION WITH ‘ INPUT EQUIVALENCES’ 


In this case, outputs are specified to certain elements ot J and, in addition, 
certain input sequences are specified as equivalent or indistinguishable (in 
the sense described in § 3), and certain other input sequences are forbidden. 
There may or may not be ‘don’t-care’ sequences, besides. 


It is to be pointed out here that if sequences i and k are equivalent, then 
for all j, O< j <2" — 1, the extended sequences il; and k1; will have to be 
equivalent. 


In forming Table I, the outputs are written in wherever they are specified 
and they are marked D for the ‘don’t-care’ sequences, if any. Distinct 
internal states are now assigned, in order, starting from S, for the input 
sequence i= 0. All the input sequences specified as equivalent are assigned 
the same internal state. The forbidden sequences are assigned the internal 
state F. Thus, if sequence i has been assigned the internal state S; and if 
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j is the next sequence (in Table I) which is not equivalent to any / for / < j, 
then j is assigned the internal state S,.,.,. 





Let M be the total number of distinct internal states assigned in making 
up Table I. Table II is now formed with (M + 2) rows and 2” columns. 
The columns refer to the input states, as usual, while the rows are named to 
correspond to the internal states, without any repetition. The first row is 
named Sp, the ith row, for 2<i<M-+1 is named §S;_, and the last row 
is named F. 

















The cells in the row S, are filled in as before. The cell (Sx, Ij), for 
1<k<M is filled in as follows: Let i(k) be the smallest input sequence 
in Table I which has been assigned the internal state S;. Form the input 
sequence i(k) I; by adding I; on the extreme right of i(k). Now, fill in the 
cell (Sx, 1;), the S — ¢ pair corresponding to the sequence i(k) I; in Table I. 
In case i (kK) I; happens to be a(N + 1)-tuple, it has to be reduced to a 1-tuple, 
for 1! <N, with the help of the given end conditions before applying the 
above rule. 


All the cells in the last row named F have F as their S-entry and ¢, as 
their ¢-entry, where ¢, is the output specified for any forbidden sequence. 
(Note that all forbidden sequences are equivalent.) 


It is a straightforward process to verify that the transition table so con- 
structed is complete and consistent, in the sense that every input sequence 
i occurring in Table I can be built out of Table II starting with S, and that 
it will terminate in the same internal state S (i) as that specified for i in Table I. 


Case 3 is now reduced to case | and the rest of the reduction process will 
go through as outlined in § 6 if there are ‘ don’t-care’ sequences, or else as 
in §5. 
































































8. CONCLUDING REMARKS 


In this paper we have been mainly concerned with the synthesis of 
sequential machines given their required behaviour in terms of events and 
associated outputs. The emphasis has been on obtaining a uniform proce- 
dure, the details of which can be programmed on a computer with ease. 
Although the formalism developed here is complete as far as it goes, it should 
be pointed out that several problems remain open. Here, we shall refer to 
only two which bear directly on the contents of this paper. 























The first is concerned with the problem of specification. Although it 
follows from Kleene’s results that everything that any finite automaton does 
can be completely described in terms of events and outputs, this does not, 
in practice, seem to be of much help. Presumably, there are behaviours 














le | 
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which are intuitively meaningful but which one does not know how to express 
in terms of events and outputs. It should be useful to have a uniform pro- 
cedure for doing this, if not in every case, at least, in a suitably well-defined 
large class of cases. 


The second problem has to do with the design of complex sequential 
machines. In practice, one tries to visualise a large machine as made up of 
submachines, each functionally independent to a large extent, and thus 
reduce the given problem to one of synthesising a large number of smaller 
machines. In such cases, there may be feedback loops from the output of 
one submachine to the input of another and problems analogous to those of 
‘ stability ’ might arise. It might be useful to consider modifications or suitable 
extensions of the formalism given here to include these cases. 


That both these problems are of great practical importance is shown by 
the fact that they play a central role in the specification and design of large-scale 
digital computers. It is of interest to note that digital computers, although 
seemingly well understood from the point of view of their description as 
sequential machines, still fall outside the scope of all existing routine synthesis 
procedures, for the most part. 


9. ILLUSTRATIVE EXAMPLES 


In this final section, we give two examples illustrating the concepts and 
the methods introduced in the paper. 


Example 1.—Synthesis of a sequential machine with one input flip-flop 
and two output flip-flops and whose output is a function of utmost the three 
previous inputs. The required behaviour is specified exhaustively by meanS 
of the assignment table in Fig. 4. (The 2 input states are denoted as 0 and 1 
and the 4 output states by 00, 01, 10 and 11 respectively.) 

End condition.—If 1,1, .... Ix-2Ik-aIk is any sequence of length k, for 
k>3, then, it is equivalent to the sequence I,_,I,_,I,. 


From Fig. 4, one obtains the modified transition table as shown in Fig. 5. 
The reduction table, at the end of complete reduction, is shown in Fig. 6. 
Thus, the equivalent states are: 

a: (So); 

a, :(S,); 

ag: (Sg, Se, Ss, Sio> Siz» S14); 

B: (Sa); 

B,: (Sg, S,, Si) 5 

y: (S55 Sp, Sis). 
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00 
00 
01 
01 
00 
10 
00 
01 
00 
10 
00 
01 
00 
10 
00 





Fic. 5. The Modified Transition Table. 


In Fig. 7 is given the state diagram of the completely reduced machine. It is 
easily seen that this diagram is not further reducible. 


Example 2.—We give here, in terms of our notation, the specification 
for the ‘ reversible counter’ considered by Huffman [pp. 276 of (2)]. 
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S; 





Sot 
S,% 
S,f 


$3", Sy2% 

S,”, 
So’, S37 

S7, 
S10”; Sia” 


a, 
Se ¢ 





Fic. 6. The Reduced Table. 





Fic. 7. The State Diagram of the Minimal Machin . 
4input states: 1), I,, I, I; 
4 ouput states: 9, $15 Pos Bs 
Restraints on input sequences: 


(1) For any Ij, lj = Ij = ]j. 
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(2) I, is forbidden. 
Output assignments: 
(1) I, has the output dp. 
(2) If i is any sequence with output ¢z,. 
Then iI, has output ¢,;,, (Mod. 4). 


il, has output ¢$9;,, (Mod. 4). 
il, has output ¢q, 


It is easily verified that, with the above, the assignment table can be com- 4 
pletely filled in and the reduction process carried through. 
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